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1ft.  Abstract 

A  theoretical  and  experimental  study  of  the  dynamic  load-deflection  curves  of  pave¬ 
ments  was  conducted  to  determine  the  dependence  of  the  measured  dynamic  stiffness 
values  of  a  pavement  on  the  type  of  vibrator  that  is  used  to  make  the  measurements, 
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be  nonlinear,  and  a  nonlinear  vibration  theory  of  pavements  is  developed  to  de¬ 
scribe  these  curves.  This  study  gives  a  method  of  determining  the  shear  modulus 
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stiffness  for  a  series  of  vibrator  baseplate  sizes.  This  method  may  be  of  prac¬ 
tical  value  for  nondestructively  determining  the  subsurface  structure  of  a  pavement. 
Volume  I,  "Experimental  Test  Results  and  Development  of  Evaluation 
Methodology  and  Procedure"  is  being  prepared  and  will  be  released  soon. 


PRICES  SUBJECT  TO  CHANGE 


13.  Kay  War4« 

Airfield  pavements 

Load  distribution 

Nondestructive  tests 

Pavements 

Vibration  tests 

IS.  Distribution  Sfotement 

Document  is  available  to  the  public 
through  the  National  Technical  Informa¬ 
tion  Service,  Springfield,  Va.  22151 

19.  Security  Ctessif.  (of  thie  report) 

Unclassified 

20.  Security  Clossif.  (of  this  pogo) 

Unclassified 

21*  No.  of  Pages 

/// 

33.  P..m 

Si 

Form  DOT  F  1700.7  ($-72)  Reproduction  of  completed  pogo  authorised 


# 

/ 


PREFACE 


•  This  study  was  conducted  during  the  period  April  1972-October 
1971*  by  personnel  of  the  U.  S.  Army  Engineer  Waterways  Experiment 
Station  (WES).  It  was  sponsored  by  the  Federal  Aviation  Administration 
through  part  of  Inter-Agency  Agreement  FA71-WAI-218,  "Development  of 
Airport  Pavement  Criteria. " 

The  study  was  done  under  the  general  supervision  of  Messrs.  J.  P. 
Sale  and  R.  G.  Ahlvin,  Chief  and  Assistant  Chief,  respectively,  of  the 
Soils  and  Pavements  Laboratory,  Mr.  R.  L.  Hutchinson,  Special  Assistant 
to  the  Soils  and  Pavements  Laboratory,  and  Mr.  H.  H.  Ulery,  Jr. ,  Chief, 
Pavement  Design  Division,  and  under  the  direct  supervision  of 
Messrs.  A.  H.  Joseph  and  J.  W.  Hall,  Jr.,  Chiefs,  respectively,  of 
the  Pavement  Investigations  Division  and  the  Evaluation  Branch.  The 
programming  for  this  study  was  done  in  part  by  Mr.  A.  P.  Park,  Soils 
Testing  Branch.  Significant  contributions  were  made  by  Mr.  J.  L.  Green, 
Evaluation  Branch.  This  report  was  written  by  Dr.  R.  A.  Weiss. 

BG  E.  D.  Peixotto  and  COL  G.  H.  Hilt  were  Directors  of  WES  during 
the  conduct  of  this  study  and  the  preparation  of  this  report.  The 
Technical  Director  was  Mr.  F.  R.  Brown. 


TABLE  OF  CONTENTS 


1.  INTRODUCTION .  13 

1.1  BACKGROUND .  13 

1.2  OBJECTIVES .  13 

1.3  SCOPE .  lit 

2.  LINEAR  OSCILLATOR  MODEL  OF  PAVEMENT  RESPONSE .  l6 

2.1  GENERAL  CONSIDERATIONS .  16 

2.2  HOMOGENEOUS  LINEAR  ELASTIC  HALF-SPACE .  16 

2.3  DYNAMIC  LOAD-DEFLECTION  CURVES .  17 

2.  It  DYNAMIC  STIFFNESS .  18 

3.  THE  NONLINEAR  MECHANICAL  MODEL .  21 

3.1  INTRODUCTION  TO  THE  NONLINEAR  MODEL .  21 

3.2  NONLINEAR  PAVEMENT-RESTORING  FORCE .  22 

3.3  EQUATION  OF  MOTION  FOR  THE  NONT, INEAR  SYSTEM .  23 

3. It  EFFECTIVE  SPRING  CONSTANT .  2k 

3.5  CALCULATION  OF  THE  DYNAMIC  STIFFNESS  AND  THE  DEFLECTION 

AMPLITUDE  FOR  NONLINEAR  PAVEMENTS .  26 

3.6  PHYSICAL  ORIGIN  OF  THE  NONLINEAR  PAVEMENT  PARAMETERS  ...  38 

3.7  CALCULATION  OF  THE  NONLINEAR  PARAMETERS .  k9 

3.8  EFFECTS  OF  THE  MECHANICAL  CHARACTERISTICS  OF  THE  VIBRATOR 

ON  THE  MEASURED  VALUES  OF  DYNAMIC  STIFFNESS .  66 

k.  EVALUATION  OF  THEORETICAL  RESULTS  OF  THE  DYNAMIC  STIFFNESS 

STUDY .  70 

4.1  EXPERIMENTAL  PROGRAM .  70 

It. 2  COMPARISON  OF  THEORETICAL  AND  EXPERIMENTAL  RESULTS  FOR 

THE  DYNAMIC  STIFFNESS .  91 

5.  SUMMARY,  CONCLUSIONS,  AND  RECOMMENDATIONS . 104 

5.1  SUMMARY . 104 

5.2  CONCLUSIONS . 104 

5.3  RECOMMENDATIONS . 106 

REFERENCES . 108 


Preceding  pege  blank 


LIST  OF  ILLUSTRATIONS 


Ft  (rare  _ Title _  Page 

1  Atypical  dynamic  response  of  the  linear  spring  model  ....  19 

2  Theoretical  dynamic  load-deflection  curves  and  dynamic 

stiffness  curves  predicted  by  the  nonlinear  spring  model  .  33 

3  Theoretical  static  load-deflection  curves  predicted  by 

the  nonlinear  spring  model . 36 

U  Theoretical  dependence  of  the  effective  static  spring 

constant  on  the  static  load . 38 

5  Frustum  of  cone  in  which  stress  and  strain  in  the  pavement 

are  assumed  to  be  confined . U3 

6  Layered  model  of  pavement  .  . . 1*7 

7  Theoretical  dependence  of  the  spring  constant,  k^  ,  on 

the  radius  of  the  \dbrator  baseplate . 59 

8  Theoretical  dependence  of  the  third-order  nonlinear  param¬ 
eter,  b  ,  on  the  radius  of  the  vibrator  baseplate  ....  60 

9  Theoretical  dependence  of  the  fifth-order  nonlinear  param¬ 
eter,  e  ,  on  the  radius  of  the  vibrator  baseplate  ....  6l 

10  Theoretical  dependence  of  the  dynamic  stiffness,  S  ,  on 

the  radius  of  the  vibrator  baseplate . 62 

11  Dynamic  load-deflection  curves  for  different  locations  on 
AC  pavements  showing  strong  pavements  to  be  more  linear 

than  weak  pavements . 73 

12  Dynamic  load-deflection  curves  for  different  locations  on 

PCC  pavements  showing  extreme  linearity  of  curves . jk 

13  Measured  values  of  dynamic  stiffness  versus  dynamic  load 

(15  Hz) . 75 

lU  Possible  experimental  evidence  for  the  existance  cf  a 

critical  frequency  value  of  *15  Hz  for  AC  pavement  ....  76 

15  Nearly  linear  dynamic  load-deflection  curves  obtained  at 

15  Hz  for  AC  pavement . 77 

16  Dynamic  load-deflection  curves  for  various  static  loads  on 

AC  pavement . 79 

17  Experimental  values  of  dynamic  stiffness  versus  static  load 

applied  by  the  WES  50-kip  vibrator  operating  at  15  Hz  on  AC 
pavement . 80 

18  Experimental  values  of  dynamic  stiffness  versus  dynamic 

load  for  various  values  of  the  static  load . 8l 


h 


1 


Figure  _ Title _ Page 


19 

Experimental  load-deflection  curves  for  a  series  of  base¬ 
plate  sizes  on  AC  pavement  . 

84 

20 

Experimented  load-deflection  curves  for  a  series  of  base¬ 
plate  sizes  on  rigid  pavement  . . 

85 

21 

Experimental  vedues  of  the  dynamic  stiffness  versus 
vibrator  baseplate  radius . . . 

the 

•  •  •  • 

86 

22 

Experimentally  derived  values  of  the  parameter 
baseplate  radius  . 

• 

o 
•  o 

versus 

ti<i 

94 

23 

Experimentally  derived  values  of  the  parameter 
baseplate  radius  . 

b 

•  • 

versus 

•  •  •  • 

94 

24 

Experimentally  derived  values  of  the  parameter 
baseplate  radius  . 

e 

•  • 

versus 

•  •  *  i 

95 

25 

Experimentally  derived  values  of  the  parameter 
vibrator  baseplate  radius . 

*0 
•  • 

versus 

•  «  •  • 

97 

26 

Experimentally  derived  values  of  the  parameter 
baseplate  radius  . 

l2 
•  • 

versus 

•  •  •  • 

97 

27 

Experimentally  derived  values  of  the  parameter 
baseplate  radius  . 

*4 
•  • 

versus 

•  •  •  • 

98 

5 


LIST  OF  TABLES 


Table  _ Title _  Page 

1  Mechanical  Characteristics  of  Vibrators . 71 

2  Comparison  of  Dynamic  Stiffness  Values  as  Measured  by  the 

Dynaflect  Vibrator  and  the  WES  l6-Kip  Vibrator  (Concrete 
Pavement ) . 87 


3  Comparison  of  Dynamic  Stiffness  Values  Measured  by  the  Road 

Rater  and  by  the  WES  l6-Kip  Vibrator  (Flexible  Pavements)  .  89 

U  Comparison  of  Typical  Dynamic  Stiffness  Values  as  Measured 

by  the  CERF  Vibrator  and  by  the  WES  l6-Kip  Vibrator  ....  90 

5  Approximate  Numerical  Values  of  the  Parameters  Appearing  in 

the  Nonlinear  Pavement  Model  (WES  l6-Kip  Vibrator) . 98 

6  Experimental  Values  of  the  Parameters  Describing  the  Sub¬ 
surface  Structure  of  Pavements  at  the  WES  Te3t  Area  ....  99 


6 


LIST  OF  SYMBOLS 


acl*ac2' 


ci 

A 


b 

B 


e 

i(wt-A) 


& 

rD(t) 

W“> 

Fp(x) 


SC 

Fv(t) 


g 

G 

■» 

J 

s 

,0 


°i*W 


.G. 

l 


h-^ »h2*^3*  *  *  * 


Contact  radius  of  vibrator  baseplate 

Critical  contact  radius  of  vibrator  baseplate 

Amplitude  of  dynamic  deflection  of  pavement  surface 
directly  beneath  vibrator  baseplate 

Third-order  nonlinear  elastic  parameter 

Function  of  Poisson's  ratio  for  the  homogeneous 
half-space 

Values  of  B  for  each  pavement  layer 

Damping  constant  of  the  pavement-vibrator  system 

Damping  constant  for  homogeneous  linear  elastic  half¬ 
space 

Algebraic  discriminant 

Fifth-order  nonlinear  elastic  parameter 

Complex  number  notation  for  a  sinusoidal  time 
dependence 

Elastic  strain  energy  density 
Dynamic  load  of  vibrator 

Magnitude  of  the  sinusoidal  dynamic  force  applied  to 
the  pavement  surface 

Pavement-restoring  force 

Static  load  of  vibrator 

Critical  static  load  of  vibrator 

Total  force  applied  to  the  pavement  surface  (static 
plus  dynamic) 

Acceleration  due  to  gravity 

Shear  modulus  of  homogeneous  elastic  half-space 
Shear  modulus  of  subgrade 

Shear  modulus  at  the  very  top  of  the  subgrade 
Shear  moduli  of  pavement  layers 
Thickness  of  the  pavement  layers 
Sums  of  pavement  layer  thicknesses 
/=! 

Effective  spring  constant  of  a  nonlinear  pavement 


7 


00 

1 


*0**2**U* * * **J 


m 

n 

Q 

*^3’ • • • 


Q 


s 


_8_ 

QG 

Q  G 
8  8 


u2,uu,u6,...u 


t 

u 

i 

V 


W 


W2*VW6 


x 

X 


*  a2  *  *  *  * 


Spring  constant  of  a  homogeneous  linear  elastic 
half-space 

Effective  quasi-static  spring  constant 

Linear  elastic  parameter  of  a  nonlinear  pavement 

Finite  depth  of  influence  of  the  static  strain  field 

Coefficients  of  the  power  series  expansion  of  the 
finite  depth  of  influence 

Lumped  mass  of  pavement  and  soil 

Number  of  layers  of  pavement 

Function  of  Poisson's  ratio 

Value  of  Q  for  each  pavement  layer 

Value  of  Q  for  the  subgrade 

Value  of  Q  at  very  top  of  subgrade 

Value  of  QG  averaged  over  the  depth  of  influence  SL 

Average  value  of  Q  G  in  the  interval  l  -  H_ 

s  s  3 

Dynamic  stiffness  of  pavement 

Critical  value  of  dynamic  stiffness 

Value  of  the  dynamic  stiffness  obtained  from  S  by 
taking  k  =  kQ  (or  A  =  0). 

Time 

Elastic  strain  energy 

Coefficients  of  the  power  series  expansion  of  the 
elastic  strain  energy 

Volume  of  the  frustum  of  a  cone  having  a  depth  equal  to 
the  finite  depth  of  influence  in. the  pavement 

Work  done  during  elastic  deflection  of  the  pavement 
under  the  action  of  the  static  load  F 

s 

Coefficients  in  the  series  expansion  of  W 

Total  elastic  deflection  of  the  pavement  surface  under 
the  vibrator  baseplate 

Velocity  of  pavement  surface 

Acceleration  of  pavement  surface 

Static  elastic  deflection  of  pavement  surface  beneath 
the  vibrator  baseplate 

Depth  below  surface  of  the  pavement 

Coefficients  appearing  in  the  power  series  expansion 
of  the  amplitude  of  the  dynamic  deflection 


8 


ei*62'**,6n 

Y 

6 

Ab 

Ae 

M 

AV 

£ 

e2»Ve6 

n 

0 

K 


k:2',cU’k6 

X 


Coefficients  appearing  in  the  power  series  expan¬ 
sion  of  the  dynamic  stiffness 

Weight  density  of  the  half-space 

Function  of  the  expansion  coefficients  of  the  fi- 
lite  depth  of  influence 

Value  of  the  discontinuity  of  b  for 

Value  of  the  discontinuity  of  e  for  £.g  = 

Increment  of  finite  depth  of  influencf 

Increment  of  the  volume  of  the  frustum  of  the  cone 
containing  the  strain  field 

Vertical  strain  in  pavement  (assumed  constant) 

Expansion  coefficients  appearing  in  the  expressions 
for  kQQ  ,  b  ,  and  e 

5/8 

Elastic  volume  dilation 

Ratio  of  the  radius  of  the  lower  base  to  the  radius 
of  the  upper  area  of  the  frustum  of  the  cone  of 
stress 

Expansion  coefficients  appearing  in  the  expressions 
for  kQ0  ,  b  ,  and  e 

Lame  elastic  constant 


u  3/b 


v  Poisson's  ratio  for  homogeneous  elastic  material 


vw 


p 

0 


02(a),03(a),01+(a) 


¥ 


Poisson's  ratio  for  subgrade 

Poisson's  ratios  for  pavement  layers 

Dynamic  elastic  deflection  of  pavement  surface 
beneath  the  vibrator  baseplate  measured  from  the 
static  equilibrium  deflection 

Function  of  the  expansion  coefficients  of  the  finite 
depth  of  influence 

Angle  which  appears  in  the  solution  of  the  cubic 
equation  for  calculating  xg  in  terms  of  Fg 

Angle  of  stress  distribution 

Functions  of  the  baseplate  radius  and  Poisson's 
ratio  of  the  successive  pavement  layers 

2  -U 

F^Sq  ,  expansion  parameter 

Volume  factor  for  the  frustum  of  the  cone 


to  Angular  frequency 


9 


Critical  angular  frequency 
Resonance  angular  frequency 

Phase  angle  between  the  dynamic  load  applied  to  the 
pavement  surface  and  the  dynamic  deflection  of  the 
pavement  surface 


10 


CONVERSION  FACTORS,  BRITISH  TO  METRIC  UNITS  OF  MEASUREMENT 


British  units  of  measurement  used  in  this  report  can  be  converted  to 
metric  units  as  follow: 


Multiply 


SL 


To  Obtain 


inches 

feet 

square  inches 

pounds  (force)  per 
inch 

pounds  (force)  per 
square  inch 

pounds  (force)  per 
cubic  inch 

kips  (force) 

kips  (force)  per 
inch 

kip-second  per 
inch 

kip-second2  per 
inch 


2.51* 

0.30U8 

6.1*516 

175.1268 


centimeters 

meters 

square  centimeters 
newtons /meter 


0.6891*757 

27,679.91 

1*, 1*1*8.222 
1.7512685 

1.7512685 

1.7512685 


newtons  per  square 
centimeters 

kilograms  per  cubic 
meter 

newtons 

kilonewtons  per 
centimeters 

kilonewtons-second  per 
centimeter 

2 

kilonewtons-second  per 
inch 


11 


1.  INTRODUCTION 


1.1  BACKGROUND 

Nondestructive  vibratory  testing  of  pavements  may  be  of  impor¬ 
tance  toward  predicting  the  performance  of  airfield  pavements  and  may 

1-3 

be  used  for  the  rapid  evaluation  of  pavement  strength.  To  be  useful, 
the  physical  quantities  measured  by  the  nondestructive  testing  technique 
must  be  related  to  pavement  performance.  Pavement  performance  is  mea¬ 
sured  by  number  of  aircraft  coverages  on  the  pavement  required  to  reach 
some  defined  condition  of  failure.  The  U.  S.  Army  Engineer  Waterways 
Experiment  Station  (WES)  was  requested  to  perform  experimental  and  theo¬ 
retical  investigations  to  determine  if  the  physical  quantities  measured 
by  the  nondestructive  technique  can  be  used  for  pavement  evaluation  and 
can  be  related  to  pavement  performance.  Some  of  the  quantities  that  are 
measurable  by  the  nondestructive  vibratory  technique  are: 

a.  The  dynamic  deflection  of  the  pavement  surface  versus  the 
frequency  of  vibratory  loading  for  a  series  of  fixed  static 
and  dynamic  loads. 

b.  The  stress  and  strain  distribution  in  the  pavement  around  the 
vibrator  measured  on  instrumented  pavement  sections. 

£.  Rayleigh  wave  dispersion  curves  giving  phase  velocity  versus 
wavelength  measured  with  the  wave  propagation  techniques. 

d.  The  dynamic  deflection  of  the  pavement  surface  versus  the 
dynamic  force  for  a  series  of  fixed  static  loads  and  fixed 
frequencies. 

Most  of  the  previous  work  on  the  nondestructive  testing  of  pavements  has 
treated  the  mechanical  quantities  listed  in  Subparagraphs  a,  b,  and  £. 
This  report  concentrates  primarily  on  the  nonlinear  response  exhibited 
by  pavements  through  the  measurements  listed  in  Subparagraph  d_.  Further 
study  of  the  physical  quantities  mentioned  in  Subparagraphs  a,  b,  and  £ 
should  be  made  in  the  light  of  the  new  results  obtained  from  the  study 
of  nonlinear  effects  in  pavements. 

1.2  OBJECTIVES 

The  overall  objectives  of  this  pavement  study  are: 

a.  The  development  of  a  mechnical  model  which  describes  the 
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measured  response  of  pavements  to  a  sinusoidal  dynamic  load¬ 
ing  that  is  applied  to  the  pavement  surface. 

b.  The  development  of  a  method  for  determining  the  subsurface 
structure  of  the  pavement  in  terms  of  the  measured  dynamic 
response  of  the  pavement. 

The  development  of  the  pavement  response  model  includes  the  fol¬ 
lowing  specific  objectives: 

a.  To  determine  the  effects  of  intrinsic  pavement  properties  and 
structure  on  the  dynamic  load-deflection  curves. 

b.  To  determine  the  effects  of  such  vibrator  characteristics  as 
dynamic  load,  static  load,  and  contact  area  on  the  dynamic 
load  deflection  curves. 

£.  To  calculate  the  dynamic  stiffness  and  determine  its  depen¬ 
dence  on  the  intrinsic  properties  of  the  pavement  and  on  the 
characteristics  of  the  vibrator  used  to  measure  this  quantity. 

d.  To  develop  a  theory  of  the  nonlinear  dynamic  response  of 
pavements  which  enables  the  comparison  of  dynamic  stiffness 
measurements  obtained  from  different  vibrators  at  the  same 
pavement  location. 

The  theoretical  work  done  in  this  report  may  have  applications 
for  the  nondestructive  testing  of  roads  and  airport  pavements.  The 
possible  practical  applications  of  this  work  are  twofold:  (a)  the  use 
of  the  dynamic  stiffness  measurement  for  determining  the  subsurface 
structure  of  the  pavement,  i.e.,  the  shear  modulus  and  thickness  for 
each  pavement  layer,  and  (b)  the  development  of  the  capability  of  com¬ 
paring  the  values  of  the  dynamic  stiffness  measured  by  different  vibra¬ 
tors  at  the  same  pavement  location. 

1.3  SCOPE 

To  achieve  the  objectives  listed  above,  both  theoretical  and  ex¬ 
perimental  studies  were  conducted. 

1.3.1  THEORETICAL  STUDIES 

The  theoretical  studies  included: 

a.  The  formulation  of  a  nonlinear  mechanical  model  to  describe 
the  response  of  a  pavement  to  static  and  dynamic  loading. 

b.  The  determination  of  effects  of  the  structure  of  the  pavement- 
soil  system  on  the  parameters  which  appear  in  the  nonlinear 
vibration  model. 


A  numerical  evaluation  of  the  parameters  that  appear  in  the 
nonlinear  model. 

d.  The  development  of  formulas  giving  the  shear  modulus  and  layer 
thickness  of  each  pavement  layer  in  terms  of  quantities  that 
are  obtained  directly  from  the  measured  dynamic  load- 
deflection  curves. 

1.3.2  EXPERIMENTAL  STUDIES 

The  experimental  studies  performed  on  both  actual  airport  pave¬ 
ments  and  especially  constructed  test  sections  included: 

a.  The  measurement  of  dynamic  load-deflection  curves  using  a 
vibrator  developed  at  WES  which  can  generate  dynamic  loads 
up  to  16  kips#  (WES  l6-kip  vibrator)  with  a  constant  l6-kip 
static  load  and  a  constant  frequency  of  15  Hz. 

b.  The  measurement  of  dynamic  load-deflection  curves  at  a  con¬ 
stant  static  load  of  l6  kips  for  a  series  of  fixed  frequen¬ 
cies  in  the  range  from  10— UO  Hz. 

£.  The  measurement  of  dynamic  load-deflection  curves  at  a  con¬ 
stant  static  load  of  16  kips  and  a  constant  frequency  of 
15  Hz  for  a  series  of  baseplates  whose  diameters  ranged  from 
5-18  in. 

d.  The  measurement  of  dynamic  load-deflection  curves  at  constant 
frequency  and  constant  baseplate  size  for  a  range  of  static 
loads  from  5-50  kips. 


#  A  table  of  factors  for  converting  U.  S.  customary  units  of  measure¬ 
ment  to  metric  (SI)  units  is  given  on  page  11. 
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2.  LINEAR  OSCILLATOR  MODEL  OF  PAVEMENT  RESPONSE 


2.1  GENERAL  CONSIDERATIONS 

Nondestructive  vibratory  testing  of  pavements  uses  a  mechanical 
vibrator  operating  at  a  known  frequency  and  dynamic  force  applied  to  the 
pavement  surface  to  produce  a  time-dependent  sinusoidal  deflection  of 
the  pavement  surface  directly  beneath  the  vibrator  baseplate.  The  mag¬ 
nitude  of  the  dynamic  deflection  of  the  pavement  surface  for  a  series 
of  dynamic  force  levels  and  frequencies  is  considered  to  be  a  measure  of 
the  strength  of  a  pavement.  This  section  discusses  a  linear  oscillator 
model  used  to  describe  the  motion  of  the  surface  of  a  linear  elastic 
half-space  and  then  shows  how  this  model  fails  to  account  for  the  mea¬ 
sured  values  of  the  dynamic  deflection  of  the  pavement  for  a  series  of 
frequencies  and  dynamic  loadings  generated  by  the  vibrator.  The  con¬ 
cepts  of  dynamic  stiffness  and  deflection  are  introduced,  and  the  sep¬ 
aration  of  static  and  dynamic  displacements  is  demonstrated. 

2.2  HOMOGENEOUS  LINEAR  ELASTIC 
HALF-SPACE 

The  equation  of  motion  of  a  mass  of  pavement  or  soil  undergoing 
vertical  oscillations  on  the  surface  of  a  homogeneous  elastic  half-space 
is 


mX  +  CHx  +  kHx  =  Fy(t) 


(2.1) 


where 


m  =  lumped  mass  of  pavement  and  soil 

X  =  acceleration  of  pavement  surface 

=  damping  constant  (3  *  4a2/Gy/g/(l  -  v),  where  a  is  the 
contact  radius  of  the  vibrator  baseplate,  G  is  the 
shear  modulus  of  the  half-space,  y  is  the  density  by 
weight  of  the  half-space,  g  is  the  acceleration  due 
to  gravity,  and  v  is  Poisson's  ratio.  Reference  4) 

x  =  velocity  of  pavement  surface 
k„  =  spring  constant  (4Ga/(l  -  v),  Reference  4) 

n 

x  =  total  elastic  deflection  of  the  pavement  surface  under 
the  vibrator  baseplate 
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Fy(t)  ■  total  force  applied  to  the  pavement  surface  (static  plus 
dynamic ) 

t  a  time 

The  values  of  k„  and  C„  are  chosen  to  construct  a  damped  spring  model 

n  it  . 

for  the  vertical  vibrations  of  an  elastic  half-space;  therefore  C„ 

n 

represents  the  radiation  damping  of  the  system.  If  viscous  friction  is 

present,  the  value  of  the  actual  damping  constant  may  be  considerably 

larger  than  the  value  of  C„  . 

xl 

The  total  force  applied  by  the  vibrator  is  written  as 

Fy(t)  =  Fg  +  FD(t)  (2.2) 

where  F  equals  static  load  and  Fn(t)  equals  dynamic  load.  The 
s  L) 

total  displacement  can  be  written  as 

x  =  xe  +  £  (2.3) 

where  xg  is  the  static  elastic  deflection  of  the  pavement  surface  be¬ 
neath  the  vibrator  baseplate,  Fg/kjj  ,  and  £  is  the  dynamic  elastic 
deflection  of  pavement  surface  beneath  the  vibrator  baseplate  measured 
from  the  static  equilibrium  deflection.  Combining  Equations  2.1,  2.2, 
and  2.3  gives  the  equation  of  motion  as 

mi  +  CRC  +  kR£  =  FD(t)  (2.1*) 

wherein  all  static  forces  and  displacements  have  canceled.  Therefore, 
for  a  linear  system,  the  static  deflection  does  not  affect  the  dynamic 
response  of  the  vibrator  mass;  only  the  reference  point  is  changed  by 
the  static  load. 

2.3  DYNAMIC  LOAD-DEFLECTION  CURVES 

For  a  sinusoidal  driving  force,  the  dynamic  deflection  obtained 
from  Equation  2.1*  is^ 
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(2.5) 


FD(w)e 


i(wt-A) 


-  ”"2) 


2  2 
+  CHW 


where : 

Fd(w)  *  magnitude  of  the  sinusoidal  dynamic  force  applied  to 
pavement  surface 

ei(wt-A)  m  coap2_ex  number  notation  for  a  sinusoidal  time  depen¬ 
dence  where  i  =  /l  ,  u  =  angular  frequency,  and  A  = 
phase  angle  between  the  dynamic  load  applied  to  the 
pavement  surface  and  the  dynamic  deflection  of  the 
pavement  surface. 

Two  kinds  of  dynamic  response  curves  of  physical  interest  can  be  ob¬ 
tained  from  Equation  2.5 : 

a.  Dynamic  deflection  versus  frequency. 

b.  Dynamic  deflection  versus  dynamic  force. 

For  a  linear  system,  the  magnitude  of  the  maximum  dynamic  deflection  is 
a  simple  linear  function  of  FD(w)  as  shown  in  Figure  la.  The  magni¬ 
tude  of  the  peak  dynamic  deflection  as  a  function  of  frequency  appears 
in  Figure  lb  for  a  constant  force  vibrator  and  for  an  eccentric  mass 
vibrator  (where  the  dynamic  force  is  frequency-dependent  in  the  manner 
Fjj(u)  ~  a)  ).  The  WES  l6-kip  vibrator  is  a  constant  force  vibrator. 
Therefore,  for  a  linear  system  the  dependence  of  £  on  u>  is  rather 
complicated,  but  the  dependence  of  £  on  Fp(u>)  is  given  simply  by  a 
straight  line  whose  slope  is  the  dynamic  stiffness.  The  phase  angle  A 
is  assumed  to  be  the  same  for  all  of  the  elements  of  the  mass  of  pave¬ 
ment  and  subgrade  which  enter  into  motion  with  the  vibrator  mass.  This 
is  the  lumped  mass  assumption,  which  requires  that  m  be  interpreted  as 
an  effective  mass  which  vibrates  in  phase  with  the  vibrator  mass  ana  has 
a  value  which  is  determined  by  requiring  that  the  theoretical  frequency 
response  curves  agree  with  the  measured  frequency  response  curves. 

2.U  DYNAMIC  STIFFNESS 

Equation  2.5  shows  that  for  a  linear  system,  the  dynamic  stiff¬ 
ness  is  given  by 
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DYNAMIC  DEFLECTION,  IN 


•.  LOAD-DEFLECTION  CURVE  FOR  LINEAR  SYSTEM 


0  5  10  15  20  25  30  35 

FREQUENCY,  HZ 

b.  FREQUENCY  RESPONSE  OF  LINEAR  SYSTEM 


Figure  1.  Typical  dynamic  response  of  the  linear  spring  model 
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(2.6) 


and  depends  only  on  the  following  quantities: 

a.  Frequency. 

b .  Spring  constant . 

c_.  Damping  constant. 

d.  Lumped  mass  of  pavement  and  soil. 

The  elastic  parameters  of  the  pavement,  G  and  v  ,  and  the  radius  of 
the  contact  area  of  the  vibrator  with  the  pavement  enter  the  dynamic 
stiffness  through  kH  and  CH  as  seen  from  Equation  2.6  and  the  ex¬ 
pressions  for  the  spring  constant  and  damping  constant.  For  a  linear 
oscillator  model,  the  dynamic  stiffness  does  not  depend  on  the  dynamic 
load  or  on  the  equilibrium  elastic  deflection,  i.e.,  £  is  a  linear 

function  of  Fp(u>)  .  However,  the  experimental  values  of  the  dynamic 
stiffness  of  pavements,  as  given  in  Section  L .  1 ,  indicate  a  strong  de¬ 
pendence  of  the  dynamic  stiffness  on  the  dynamic  load  and  on  the  static 
elastic  equilibrium  displacement  of  the  pavement  surface.  Therefore, 
the  linear  oscillator  model  is  insufficient  to  describe  the  response  of 
pavements  to  dynamic  loadings,  and  a  nonlinear  oscillator  model  is  re¬ 
quired  to  explain  the  experimental,  data. 
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3.  THE  NONLINEAR  MECHANICAL  MODEL 


3.1  INTRODUCTION  TO  THE  NONLINEAR 
MODEL 

In  this  section,  a  nonlinear  mechanical  model  is  developed  to 
describe  the  response  of  a  pavement -subgrade  system  to  a  sinusoidal  dy¬ 
namic  loading  applied  to  the  surface  of  the  pavement  by  a  vibrator.  The 
model  is  developed  in  three  basic  steps: 

a.  The  determination  of  the  nonlinear  pavement -restoring  force 
in  terms  of  three  parameters:  the  linear  elastic  parameter 
of  a  nonlinear  pavement,  the  third-order  nonlinear  elastic 
parameter,  and  the  fifth-order  nonlinear  elastic  parameter. 

b.  The  solution  of  the  motion  equation  (2.1*)  for  the  case  of  the 
nonlinear  pavement-restoring  force  and  the  subsequent  cal¬ 
culation.  of  the  dynamic  stiffness  and  deflection  of  the  pave¬ 
ment  as  a  xVinction  of  the  static  and  dynamic  loads  exerted  by 
the  vibrator. 

c_.  The  determination  of  the  parameters  Icqq  ,  b  ,  and  e  in 
terms  of  the  elastic  constants  of  the  layered  pavement- 
subgrade  system  and  in  terms  of  the  finite  depth  of  influence 
that  a  static  surface  load  produces  in  this  system. 

If  for  a  fixed  frequency  the  dynamic  deflection  of  the  pavement 
surface  is  not  directly  proportional  to  the  dynamic  force,  the  system  is 
said  to  be  nonlinear.  The  experimental  data  of  Section  U .  1  indicate 
that  this  is  the  case  for  most  asphaltic  concrete  (AC)  pavements  and  for 
some  Portland  cement  concrete  (PCC)  pavements.  It  will  be  shown  in 
Section  3-8.2  that  the  nonlinear  behavior  of  a  pavement  undergoing 
forced  sinusoidal  vibrations  can  produce  very  different  values  of  dy¬ 
namic  stiffness  such  as  those  measured  at  the  same  location  by  differ¬ 
ent  mechanical  vibrators.  Therefore,  it  is  important  to  be  able  to 
account  for  the  nonlinear  effects  by  a  simple  physical  model. 

A  physical  and  mathematical  model  for  the  nonlinear  response  of 
pavements  can  be  derived  which  will  account  for  the  dependence  of  the 
impedance  values  on  the  type  of  vibrator  used  to  determine  them,  i.e., 
on  the  static  weight,  dynamic  load,  and  contact  area  of  the  vibrator. 
This  report  will  show  that  it  is  possible  to  describe  the  dependence  of 
the  measured  values  of  pavement  dynamic  stiffness  on  the  physical 
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characteristics  of  the  vibrator  by  introducing  three  parameters  to  de¬ 
scribe  the  nonlinear  pavement-restoring  force. 

3.2  NONLINEAR  PAVEMENT-RESTORING 
FORCE 

The  pavement-restoring  force  is  the  force  that  the  bulk  pavement 
exerts  on  the  lumped  pavement  mass  from  below.  In  linear  Equation  2.1, 
the  pavement-restoring  force  is  simply  k^x  .  In  general ,  the  pavement¬ 
restoring  force  is  not  equal  to  the  force  generated  by  the  vibrator; 
only  for  the  static  case  are  these  two  forces  equal.  The  first  task 
to  be  accomplished  is  the  development  of  a  mathematical  expression  for 
the  pavement-restoring  force  which  satisfies  the  following  two  very 
general  criteria: 

a.  The  mathematical  form  of  the  pavement -re storing  force  will 
be  sufficiently  general  so  that  the  nonlinear  dynamic  re¬ 
sponse  of  the  pavement  that  is  calculated  from  this  restoring 
force  will  be  adequate  to  describe  the  experimental  nonlinear 
dynamic  load-deflection  curves. 

b.  Only  terms  based  on  sound  physical  theory  are  included  in 
the  mathematical  form  of  the  pavement-restoring  force. 

The  form  of  the  nonlinear  elastic  pavement-restoring  force  used 
in  the  nonlinear  model  is  determined  by  requiring  the  restoring  force 
to  be  antisymmetric  in  the  deflection  of  the  pavement  surface,  i.e., 

Fp(x)  =  -Fp(-x)  (3.1) 

where  Fp(x)  equals  the  pavement-restoring  force.  Equation  3.1  is  sat¬ 
isfied  for  the  linear  case,  Fp  = 
restoring  force  which  satisfies  Equation  3.1,  and  which  is  found  to  be 
adequate  to  describe  the  dynamic  load-deflection  curves  for  pavements,  is 

Fp(x)  =  kQOx  +  bx3  +  ex5  (3.2) 

where  k^  equals  the  linear  elastic  parameter  of  a  nonlinear  pavement 
while  b  and  e  equal  respectively  the  third-  and  fifth-order  non¬ 
linear  parameters.  The  experimental  data  of  Section  Ij.1  indicate  that 
at  least  two  nonlinear  elastic  parameters,  b  and  e  ,  are  required  to 
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A  simple  nonlinear  pavement- 
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describe  AC  and  PCC  pavements.  The  linear  ;;pring  constant  kQ0  which 
appears  in  Equation  3.2  is  not  in  general  equal  to  the  spring  constant 
k„  which  describes  the  homogeneous  linear  elastic  half-space. 

n 

3.3  EQUATION  OF  MOTION  FOR  THE 
NONLINEAR  SYSTEM 

The  equation  of  motion  of  the  oscillating  lumped  pavement  and 
soil  mass  can  now  be  written  using  the  expression  for  the  nonlinear 
pavement-restoring  force  derived  in  the  previous  section.  This  equation 
of  motion  cannot  be  completely  separated  into  static  and  dynamic  parts 
as  was  the  case  for  the  linear  elastic  system. 

The  equation  of  motion  for  the  nonlinear  spring  is  given  by 

mX  +  Cx  +  k0Qx  +  bx3  +  ex5  =  F  (t)  (3.3) 

where  C  is  the  damping  constant  of  the  pavement-vibrator  system,  and 
m  is  the  in-phase  lumped  mass  of  the  pavement  and  subgrade.  The  value 
of  C  is  larger  than  the  value  of  the  radiation  damping  constant  CH 
which  appears  in  Equation  2.6  because  C  describes  several  material 
damping  processes  in  addition  to  the  dissipation  of  energy  by  mechanical 
radiation.  Equation  3.3  can  be  greatly  simplified  by  choosing  a  new 
origin  of  coordinates  as  in  Equation  2.3,  such  that  the  motion  is  de¬ 
scribed  in  terms  of  coordinates  measured  relative  to  the  static  equilib¬ 
rium  deflection.  By  itself  the  static  load  produces  a  static  deflection 
given  by 

F  =  kAAx  +  bx3  +  ex5  (3.10  ' 

s  00  e  e  e 

Substituting  Equation  3.1*  into  Equation  3.3  enables  the  equation  of 
motion  to  be  written  as 

mX  +  Cx  +  kQ0(x  -  xe)  +  b(x3  -  x3)  +  e(x5  -  x^)  =  F])(t)  (3.5) 

Using  Equation  2.3  and  the  fallowing  algebraic  identities: 

x3  -  x3  =  (x  -  x  ) (x^  +  xx  +  x^)  (3.6) 

e  e  \  e  e  / 
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x5  -  x^  ■  (x  -  xfi)(x1+  +  x3xr  +  x2x2  +  XX3  +  xr)  (3.7) 


allows  Equation  3.5  to  be  rewritten  as 


m6  +  Cl  +  kQ5  +  b?3  +  ef-5  =  FD(t)  (3.8) 

where  kQ  is  the  effective  quasi-static  spring  constant  and  is  defined 
by 


ko  s  koo 


+  3bx2  +  5exR  +  g(xRS) 


(3.9) 


and 


g(x  e)  =  3bx  £  +  10ex3?  +  10ex2£2  +  5ex  £3  (3.10) 

e  e  e  e  e 

Equation  3.8  is  a  generalization  of  the  Duffing  Equation.^ 

3.U  EFFECTIVE  SPRING  CONSTANT 

In  this  section,  the  equation  which  determines  the  amplitude  of 
the  sinusoidal  dynamic  deflection  of  the  pavement  surface  beneath  the 
vibrator  mass  is  developed.  The  amplitude  equation  is  expressed  in 
terms  of  an  effective  spring  constant  which  in  turn  depends  on  the 
static  and  dynamic  deflections  of  the  pavement  surface.  The  dynamic 
stiffness  for  the  nonlinear  system  will  eventually  be  expressed  in  terms 
of  this  effective  spring  constant. 

The  functions  kQ(xeC)  and  g(xg£)  are  time-dependent,  and 
therefore  Equation  3.8  is  very  difficult  to  solve  exactly.  Under  spe¬ 
cial  conditions  to  be  described,  the  coefficient  which  appears  in  Equa¬ 
tion  3.8  may  be  taken  to  be  independent  of  time,  thereby  making  this 
equation  somewhat  easier  to  solve.  For  harmonic  motion,  the  dynamic 
force  applied  to  the  pavement  surface  by  the  vibrator  can  be  written  as 

FD(t)  =  FD(u>)eia)t  (3.11) 


2k 
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£(t)  »  Aei(ut‘A)  (3.12) 

where  A  equals  the  amplitude  of  the  dynamic  deflection  of  the  pavement 
surface  directly  beneath  the  vibrator  baseplate.  The  dynamic  deflection 
of  the  lumped  mass  is  assumed  to  be  equal  to  the  dynamic  deflection  of 
the  pavement  surface.  For  the  case  in  which  the  dynamic  deflection 
amplitude  is  much  less  than  the  static  equilibrium  deflection,  A  «  xg  , 
g(x  £)  **  0  can  be  used  in  Equation  3.9,  while  for  the  case  where 
A  **  xg  ,  the  time-averaged  value  of  g(xe£)  **  10ex|[  can  u8e<*  *n 
Equation  3.7.  For  the  two  special  cases,  the  coefficient  kQ  can  be 
written  as 

k0  =  k00  +  3bxe  +  5exi  A  <<  xe  (3.13) 

k0  =  k00  +  3bx5  +  15eXe  A*Xe  {3’lk) 

A  simple  linear  interpolation  formula  for  kQ  is  given  by 

k0  =  k00+  3(1+  <3'15> 

It  should  be  noted  that  the  choice  of  k^  as  time-independent  is  an  ap¬ 
proximation  which  becomes  invalid  for  large  dynamic  deflections. 

Even  with  the  coefficient  k^  assumed  to  be  time-independent. 
Equation  3*8  is  a  nonlinear  equation.  However,  it  can  be  shown  that 
Equation  3.8  can  be  cast  into  the  form  of  an  equivalent  linear  system 
whose  amplitude  equation  is  ^ 

A2[(k  -  mu2)  +  C2w2]  =  F2(u>)  ( 3 •  16 ) 

provided  an  effective  spring  constant  is  introduced  which  is  defined  by 

k  =  kQ  +  pbA2  +  neAU  (3.17) 

where  k^  is  given  by  Equation  3.15,  P  =  3/b  ,  and  n  =  5/8  .  The 
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effective  spring  constant,  k  ,  is  seen  to  be  a  function  of  the  ampli¬ 
tude  of  the  dynamic  deflection  and  also  depends  on  the  static  equilib¬ 
rium  deflection  through  the  coefficient  kQ  .  If  the  static  load  ap¬ 
plied  to  the  pavement  by  the  vibrator  were  zero,  then  xg  *  0  , 
g(xg5)  *  0  ,  and  k^  *  k^  .  For  this  case,  there  could  be  no  coupling 
of  terms  between  xfi  and  £  and  the  effective  spring  constant  would 
depend  only  on  the  amplitude  of  the  dynamic  deflection.  On  the  other 
hand,  if  the  dynamic  load  were  zero,  the  effective  spring  constant  would 
be  k  =  kQ  and  would  depend  only  on  the  static  equilibrium  deflection. 

3.5  CALCULATION  OF  THE  DYNAMIC  STIFF¬ 
NESS  AND  THE  DEFLECTION  AMPLI¬ 
TUDE  FOR  NONLINEAR  PAVEMENTS 

This  section  considers  the  calculation  of  the  dynamic  stiffness 
and  the  dynamic  deflection  of  the  pavement  surface  and  dynamic  forces 
generated  by  the  vibrator.  The  deflection  amplitude  equation  (3.16) 
derived  in  the  previous  section  is  expanded  in  powers  of  the  deflection 
amplitude,  A  ,  to  give  a  tenth-order  algebraic  equation  for  the  de¬ 
termination  of  A  .  Infinite  series  expansions  for  the  dynamic  ampli¬ 
tude  and  the  dynamic  stiffness  are  obtained  as  solutions  to  this  equa¬ 
tion.  These  solutions  express  the  dynamic  stiffness  and  deflection  as 
functions  of  the  dynamic  load  generated  by  the  vibrator  and  the  static 
deflection  of  the  pavement  surface.  The  static  deflection  is  then  ex¬ 
pressed  in  terms  of  the  static  load,  so  that  finally  the  dynamic  stiff¬ 
ness  and  deflection  are  expressed  in  terms  of  the  static  and  dynamic 
loads  at  which  the  vibrator  is  operated. 

3-5.1  EQUATION  FOR  THE  DYNAMIC 
DEFLECTION  AMPLITUDE 


The  explicit  equation  for  the  dynamic  deflection  amplitude  will 
now  be  calculated.  The  dynamic  stiffness  of  the  pavement  which  is  de¬ 
scribed  by  a  nonlinear  oscillator  is  given  by 


S 


2 


P2  2 
C  a) 


(3.18) 


where  the  effective  spring  constant  is  given  by  Equation  3. IT.  The 
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dynamic  stiffness  depends  on  the  amplitude  of  the  dynamic  displacement 
and  the  static  equilibrium  deflection.  The  amplitude  of  the  dynamic 
deflection  is  determined  by  Equation  3.16,  which  may  be  written  as 


A2S2  ■  Fp(w) 


(3.19) 


Using  Equations  3.17  and  3.18,  the  amplitude  equation  (3.19)  can  be 
written  as 


SqA2  +  2(Jb^kQ  -  *■  ^2n(kg  -  mw2)e  +  w^b”1  J  A^ 


+  2wnbeA8  +  n2e2A10  =  F2(u>)  (3.20) 


where  SQ  is  the  value  of  the  dynamic  stiffness  obtained  from  S  by 
taking  k  =  kQ  (or  equivalently  A  =  0  in  Equation  3.17)  and  is 
defined  by  the  equation: 


S 


2 

0 


t  p2  2 
+  C  a) 


(3.21) 


Whereas  the  simple  linear  system  produces  a  linear  equation  for  the 
calculation  of  the  dynamic  displacement  in  terras  of  the  dynamic  force, 
the  nonlinear  system  appropriate  to  describe  dynamic  pavement  response 
produces  a  fifth-order  equation  for  calculating  A  in  terms  of  Fp(u>)  • 
The  value  appearing  in  Equation  3.21  is  the  dynamic  stiffness  in 

the  limit  of  zero  dynamic  loading. 

3.5.2  POWER  SERIES  EXPANSION 

The  tenth-order  equation  (3.20)  will  now  be  solved  for  the  dy¬ 
namic  amplitude  A  which  will  take  the  form  of  an  infinite  series  ex¬ 
pansion.  The  dynamic  stiffness  is  calculated  in  terms  of  A  by  Equa¬ 
tion  3.19  30  that  S  also  will  have  the  form  of  an  infinite  series 
expansion. 

The  solution  of  Equation  3.20  for  the  amplitude  of  the  dynamic 
displacement  in  terms  of  the  amplitude  of  the  dynamic  force  is,  in 
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general,  difficult  to  obtain  analytically.  For  the  case  in  which  the 
dynamic  force  is  not  very  large,  the  amplitude  of  motion  and  the  dynamic 
stiffness  are  easily  obtained  from  Equation  3.20  in  the  form  of 


Vw)  /  2  \ 

A  ■  g-  \1  +  +  a2i|i  +  . . .  j 

0(i  +  Sj*.  s/  ♦ 


S  ■  S 


(3.22) 

(3.23) 


where 

a^,  Og  =  coefficient  appearing  in  the  power  series  expansion  of 
the  amplitude  of  the  dynamic  deflection 

\p  -  expansion  parameter 

@1*  =  coefficients  appearing  in  the  power  series  expansion  of 

the  dynamic  stiffness 


The  values  of  ip  ,  ,  and  6^  can  be  obtained  by  com¬ 

bining  Equations  3.20  and  3.22  with  t.ie  following  results: 


*  = 


fd<“> 


(3.2M 


and 


a 


2 


(3.25) 

(3.26) 

(3.27) 

(3.28) 


The  solutions  given  in  Equations  3.22  and  3.23  are  valid  provided  the 
dynamic  load  is  not  so  large  as  to  prevent  the  convergence  of  these 
power  series  solutions.  Equations  3.22-3.23  have  been  derived  from 
Equation  3.19  and  give  the  fundamental  description  of  the  nonlinear  dy¬ 
namic  load-deflection  curves.  These  equations  will  be  fitted  to 
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experimental  dynamic  load-deflection  curves  in  Section  k.2.2. 

Equation  3.22  shows  that  the  amplitude  of  the  dynamic  deflection 
is  not  a  linear  function  of  but  approaches  the  linear  condition  for 

small  Fd  or  large  SQ  .  The  linear  system  can  be  regained  by  setting 
b  =  e  =  0  in  Equations  3.22-3.28.  When  =  0  ,  A  =  0  and  S  *  , 

and  the  dynamic  stiffness  depends  only  on  the  static  equilibrium  deflec¬ 
tion.  If  F  =  0  ,  the  deflection  and  stiffness  are  given  by  Equations 
s 

3.22-3.28  with  the  provision  that  kQ  be  replaced  by  the  constant  kQQ  . 
The  static  equilibrium  deflection  can  always  be  expressed  in  terms  of 
the  static  load  through  Equation  3.1*.  Therefore,  in  general,  the  dy¬ 
namic  stiffness  of  a  nonlinear  system  will  depend  on  the  magnitude  of 

F  and  F  .  The  dependence  of  S  on  F  enters  through  the  expan- 
US  u 

sion  parameter  ip  given  in  Equation  3.2i*,  while  the  dependence  of  S 

on  F  enters  through  the  function  Sn  given  by  Equation  3.21. 
s  u 

3.5.3  FIRST-ORDER  NONLINEAR  TERM 

The  expression  for  the  dynamic  amplitude  A  given  in  Equation 

3.22  shows  that  A  does  not  depend  linearly  on  F  .  The  departure 

2  D 

from  linearity  is  due  to  the  terms  ,  a \p  ,...,  that  appear  in 

Equation  3.22.  It  is  desirable  to  determine  the  physical  quantities 
which  determine  the  degree  of  departure  from  the  linear  condition, 

A  =  Fd/Sq  .  In  the  range  of  small  F^  ,  the  predominant  term  describing 
the  nonlinear  behavior  of  the  deflection  of  the  pavement  mass  is  ob¬ 
tained  from  Equations  3.22,  3.2H,  and  3.25  as  follows: 


In  general  the  degree  of  nonlinearity  depends  on  four  quantities: 

a.  The  magnitude  of  the  nonlinear  parameters  b  and  e  . 

b.  The  relative  magnitudes  of  F  and  F  . 

"  s  u 

c_.  The  frequency  at  which  the  vibrator  is  operated. 

<1.  The  static  stiffness  Sq  of  the  pavement-vibrator  system. 
The  parameter  ip  ,  which  appears  in  the  infinite  series  expansion  for 
A  and  S  in  Equations  3.22-3.2^,  depends  inversely  on  in  the 
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-1* 

form  Sq  ,  and  therefore  it  follows  that  the  dynamic  load-deflection 
curves  of  stiff  pavements  are  more  linear  than  those  of  the  more  flexi¬ 
ble  pavements.  Thus  concrete  pavements  are  expected  to  have  a  more 
linear  response  to  a  dynamic  loading  than  do  the  more  flexible  asphalt 
pavements.  The  value  of  Sq  includes  the  effects  of  the  subgrade  as 
well  as  the  effects  of  each  layer  in  the  pavement. 

3.5-1*  CRITICAL  FREQUENCY 


It  is  clear  from  Equations  3.22  and  3.29  that  the  degree  of  de¬ 
parture  from  the  linear  condition  expected  for  the  response  of  a  pave¬ 
ment  to  an  applied  vibratory  load  at  the  pavement  surface  depends  on 
the  frequency  at  which  the  vibrator  is  operated.  In  particular,  it  is 
apparent  from  Equation  3.29  that  the  first-order  nonlinear  term  is 

frequency-dependent  and  that  this  first-order  term  will  vanish  at  a 

2 

special  critical  frequency  for  which  kg  -  mui  =  0  . 

It  is  a  characteristic  property  of  the  first-order  nonlinear 
term  (Equation  3.29)  that  there  is  a  critical  frequency  for  which  this 
term  vanishes;  the  critical  frequency  is  defined  by 


(3.30) 


where  u>c  is  the  critical  angular  frequency.  In  terms  of  the  critical 
frequency,  the  first-order  nonlinear  coefficient  can  be  written  as 


a 


1 


(3.31) 


At  the  critical  frequency,  the  departure  from  a  linear  system  occurs 

2 

only  through  the  second-order  and  higher  terms  in  ,  i.e.  ,  a  ip 

3  d 

+  m  +....  Therefore  at  the  critical  frequency,  the  pavement  response 

for  small  dynamic  loads  should  be  nearly  linear.  The  critical  frequency 
depends  on  the  vibrator  characteristics  as  well  as  on  pavement  proper¬ 
ties.  The  connection  between  the  resonance  frequency  and  the  critical 
frequency  is  obtained  from  Equation  3.17  as  follows: 
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(3.32) 
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2  +  HbA2  +  _  2/C_\2 

c  m  \2m/ 


where  wn  is  the  resonance  angular  frequency.  In  general  u.,  <  u> 
n  n  c 


3.5*5  CRITICAL  STATIC  FORCE 


In  addition  to  a  critical  frequency,  uc(^a)  which  depends  on 
the  static  load  of  the  vibrator,  there  is  a  critical  static  load  F 

sc 

for  each  operating  frequency  of  the  vibrator,  which  is  defined  from 
Equation  3.29  by  =  0  or 


k.(F  )  =  mu 
0  sc7 


(3.33) 


Using  Equations  3.15,  3.31,  and  3.33,  the  first-order  coefficient 
can  be  written  as 


( 3. 3U ) 

(3.35) 


It  is  possible  to  operate  a  vibrator  at  the  critical  condition  by  ad¬ 
justing  either  the  frequency  or  the  static  load  of  the  vibrator. 


3.5-6  SHAPE  OF  THE  DYNAMIC 
LOAD-DEFLECTION  CURVES 


An  approximately  linear  dynamic  deflection  versus  dynamic  force 

curve  occurs  at  the  critical  frequency.  For  an  arbitrary  frequency,  the 

departure  from  this  approximately  linear  curve  is  positive  or  negative 

depending  on  the  sign  of  the  parameter  in  Equations  3.22  and  3.31. 

The  sign  of  the  parameter  depends  on  the  sign  of  the  parameter  b 

and  whether  w  ^  u>  or  F  ^  F  .It  will  be  shown  in  Section  3.7 
c  s  sc 

that  b  is  generally  negative.  For  the  case  b  <  0  ,  which  corresponds 
to  the  case  in  which  the  shear  modulus  of  the  half-space  is  constant 
with  depth,  or  to  the  case  of  a  layered  system  which  has  G  decreasing 
with  depth,  as  is  usually  the  case  with  pavements,  the  dynamic  stiffness 
and  deflection  versus  dynamic  force  curves  are  shown  schematically  in 
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Figure  2a.  For  the  case  b  >  0  ,  which  may  be  possible  for  the  situa¬ 
tion  in  which  the  shear  modulus  of  a  layered  system  increases  rapidly 
with  depth,  the  dynamic  stiffness  and  displacement  versus  dynamic  force 
curves  are  shown  schematically  in  Figure  2b.  Therefore  the  choice  of 
sign  b  <  0  has  physical  relevance  to  pavement  problems.  For  the 
choice  b  <  0  ,  the  sign  of  the  parameter  a  can  be  positive  or  nega¬ 
tive  depending  on  whether  ui  <  u>c  or  u  >  w  ,  respectively.  The  value 
of  u)c  can  be  determined  by  observing  the  frequency  which  produces  the 
most  linear  load-deflection  curve.  The  algebraic  signs  of  b  ana  e 
can  be  determined  from  the  manner  in  which  the  dynamic  load-deflection 
curves  bend  away  (as  in  Figures  2a  and  2b)  from  the  approximately  linear 
load-deflection  curve  which  occurs  at  u>  =  w  .  It  should  be  pointed 
out  that  the  definition  of  dynamic  stiffness  as  being  equal  to  the  ratio 
of  the  dynamic  load  to  the  dynamic  deflection  for  each  point  on  the  dy¬ 
namic  load-deflection  curve  is  different  from  that  used  in  Volume  I  of 
this  report  in  which  the  dynamic  stiffness  is  taken  to  be  a  single  num¬ 
ber  determined  from  the  portion  of  the  dynamic  load-deflection  curve 
where  F^  is  large. 

3.5-7  STATIC  EQUILIBRIUM  DISPLACEMENT 

The  explicit  dependence  of  the  dynamic  stiffness  on  is  given 

by  Equations  3.22-3.28.  These  equations  will  also  give  the  explicit 
dependence  of  S  on  Fg  ,  provided  that  the  static  equilibrium  dis¬ 
placement  x  is  expressed  explicity  as  a  function  of  F  by  using 

6  S 

Equation  3.^.  Because  Equation  3.^  is  an  equation  of  fifth  degree, 

numerical  methods  are  generally  required  for  its  solution.  However,  in 

the  extremes  of  very  large  and  very  small  values  of  F  ,  analytical 

s 

solutions  of  this  equation  are  possible.  For  a  very  small  static  load, 
the  equilibrium  elastic  displacement  is  given  by 


For  somewhat  larger  values  of  Fg  ,  the  cubic  term  manifests  itself  and 

x  may  be  obtained  from  the  approximate  equation: 
e 
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(3.37) 


3  k00  Fs  n 

x  +  — r—  x  -  —  =  0 
e  Deb 


The  discriminant  of  this  cubic  equation  is 


D  = 


kb£ 


00 

27b” 


(3.38) 


and  is  negative  for  small  F 


when  b  <  0  .  For  the  condition  D  <  0 
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the  solution  of  the  cubic  equation  can  be  written  as 

Tk 
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cos  0  =  - 


cos 
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(3.39) 

(3. bO) 


where  0  is  the  angle  which  appears  in  the  solution  of  this  cubic  equa¬ 
tion. 

Equation  3.39  has  the  value 


In  the  limit  of  small  F  (or  small  b),  the  cosine  term  in 
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(3.141) 


Combining  Equations  3.39  and  3 - Ul  gives 


x  =  — 
e  k 


00 


bF' 


00 


(3.142) 


The  solutions  of  Equations  3.39  and  3.142  have  been  derived  for  b  <  0 
and  are  therefore  applicable  to  pavements.  It  can  be  shown  that  Equa¬ 
tion  3. 1+2  is  also  valid  for  b  >  0  . 

With  increasing  values  of  F^  ,  the  fifth-order  terms  become  dom¬ 
inant,  and  in  this  region  the  approximate  solution  for  xg  is 


■&-*er-s(8T 


( 3. 1*3) 
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Equation  3.1*  is  easily  solved  for  the  general  case  of  an  arbitrary  value 
of  F  by  using  a  digital  computer.  A  schematic  graph  of  x  versus 

S  6 

F  for  pavements  (b  <  0)  is  given  in  Figure  3a  while  the  corresponding 
s 

graph  for  a  b  >  0  formation  is  given  in  Figure  3b. 


3.5.8  DEPENDENCE  OF  THE  SPRING  CONSTANT 
ON  THE  STATIC  FORCE 

The  spring  constant  k^  given  by  Equations  3.13-3.15  has  a  con¬ 
ventional  interpretation  only  when  the  dynamic  deflection  amplitude  sat¬ 
isfies  A  <<  xg  .  For  A  ^  0  the  spring  constant  has  the  approximate 
value  given  by  Equation  3.15.  Using  Equations  3.13,  3.1*2,  and  3.1*3,  the 
value  of  kg  for  zero  dynamic  amplitude  and  for  small  Fs  is 


(3.1*5) 


For  very  small  Fg  ,  Equation  3.1*1*  can  be  rewritten  as 


while  for  very  large  Fg  ,  Equation  3.1*5  can  be  rewritten  as 


(3.»*6) 


(3.1*7) 


Equations  3.1*  and  3.15  are  easily  solved  simultaneously  on  a  digital 
computer  to  give  the  general  solution,  k^  =  kg(Fg).  Figures  ha  and  b 
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show  respectively  the  dependence  of  k  on  F  for  pavements  (b  <  0) 

u  s 

and  for  b  >  0  formations.  For  the  case  b  <  0  ,  the  function  kQ 

exhibits  a  local  minimum  value  for  some  value  of  F  (or  x  ). 

s  e 

3.5.9  DAMPING  CONSTANT 

The  dynamic  stiffness  defined  by  Equation  3.18  depends  on  the 
damping  constant  C  of  the  pavement  as  well  as  on  the  effective  spring 
constant  k  .  This  damping  constant  is  not  in  general  equal  to  the 
damping  constant  for  the  homogeneous  linear  elastic  half-space,  C..  , 
that  was  defined  in  the  damping  constant  expression.  A  theoretical  cal¬ 
culation  of  C  was  not  made  in  this  report;  however,  the  nonlinear 
elastic  nature  of  flexible  pavements  gives  rise  to  a  simple  method  of 
estimating  the  value  of  the  damping  constant.  Equations  3.17,  3.18,  and 
3.30  show  that  when  u>  s  u>  ,  the  dynamic  stiffness  has  the  critical 
value : 

Sc  =  ^(ybA2  +  qeA^ )  +  C2ui2 

*  >/2Cir  (3.W) 

Therefore  a  measure  of  the  damping  constant  can  be  determined  directly 
from  the  dynamic  load  versus  deflection  curves  by  measuring  the  critical 
value  of  the  dynamic  stiffness.  An  approximation  to  the  value  of  the 
damping  constant  is  thus  given  by 


3.6  PHYSICAL  ORIGIN  OF  THE  NONLINEAR 
PAVEMENT  PARAMETERS 

3.6.1  INTRODUCTION  OF  THE  FI¬ 
NITE  DEPTH  OF  INFLUENCE 

This  part  of  the  report  examines  the  physical  origin  of  the  pa¬ 
rameters  k^  ,  b  ,  and  e  that  were  introduced  in  Section  3.2  as  part 
of  a  model  developed  to  account  for  the  observed  nonlinear  vibration 
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data  that  will  be  presented  in  Section  4.1.  It  is  the  departure  from 
perfect  linear  elasticity  that  is  responsible  for  the  nonlinear  response 
of  vibrators  operating  at  the  surface  of  pavements  and  soil  formations. 
Over  a  large  depth  in  the  half-space,  the  departure  from  a  perfectly 
linear  constitutive  equation  for  soils  and  pavements  has  the  effect  of 
introducing  a  finite  range  of  influence  of  the  elastic  stress  and  strain 
due  to  a  static  surface  loading  on  an  elastic  half-space.  In  this  re¬ 
port  the  finite  range  of  influence  of  the  static  stress  and  strain  is 
assumed  to  produce  the  experimentally  observed  nonlinear  dynamic  load 
deflection  curves. 

3.6.2  DEPARTURE  FROM  LINEAR  ELASTICITY 

Many  authors  have  studied  the  problem  of  the  static  embedding  of 

g 

a  circular  punch  into  a  linear  elastic  half-space.  Layered  elastic 

half-spaces  have  been  considered  and  so  has  the  case  in  which  G  has  a 

8“10 

power  law  variation  with  depth.  -  The  case  in  which  G  has  a  simple 
linear  variation  with  depth  has  been  worked  out  completely.11  All  of 
these  theories  are  insufficient  to  explain  the  experimental  data  appear¬ 
ing  in  Section  4.1  because  these  theories  are  based  on  the  assumption  of 
a  perfect  linear  elasticity.  One  of  the  characteristics  of  a  nonhomo- 
geneous  perfectly  linear  elastic  half-space  is  that  the  static  load- 
deflection  curve  and  the  dynamic  load-deflection  curve  are  straight 
lines.  For  instance,  a  stack  of  perfectly  linear  elastic  layers  of 
different  values  of  v  and  G  produces  a  linear  load-deflection  curve 
for  a  load  applied  to  the  surface  of  the  stack.  This  is  so  because  the 
static  stress  and  strain  in  a  perfectly  linear  elastic  half-space  ex¬ 
tend  to  an  infinite  depth  and  to  an  infinite  radial  distance,  and  the¬ 
oretically  the  entire  stack  of  elastic  layers  contributes  to  the  deflec¬ 
tion  of  the  load  at  the  surface  of  the  stack.  In  a  perfectly  linear 
theory  of  elasticity,  the  static  load-deflection  curve  for  a  stack  of 
layers  can  always  be  reproduced  by  an  equivalent  homogeneous  three- 
dimensional  half-space  with  an  equivalent  spring  constant.  Therefore, 
a  nonhomogeneous  linear  elastic  system  exhibits  a  linear  load-deflection 
curve,  and  nonhomogeneity  by  itself  is  insufficient  to  explain  the 


39 


nonlinear  load-deflection  curves  observed  on  pavements  and  soils.  A 
basic  departure  from  linear  elasticity  is  required  to  account  for  the 
observed  nonlinear  response  of  vibrators  which  are  operated  on  the  sur¬ 
face  of  airfield  pavements.  The  nonlinear  behavior  of  pavements  is  an 
example  of  the  inadequacy  of  linear  elasticity  to  explain  the  dynamic 
response  of  real  materials. 

The  static  stress  and  strain  field  has  a  finite  range  of  influ¬ 
ence  in  real  systems  such  as  PCC,  AC,  or  soil.  This  is  a  fundamental 
characteristic  of  real  media.  If  a  static  load  is  placed  on  the  surface 
of  a  soil  formation  or  a  pavement,  the  static  stress  and  strain  attenu¬ 
ate  rapidly  with  depth  and  radial  distance  and  go  to  zero  at  some  finite 
depth  below  the  surface  and  at  some  finite  radial  distance  from  the 
load.  The  stress  and  strain  field  in  a  real  half-space  does  not  extend 
to  infinite  depth  as  it  does  for  the  case  of  the  Boussinesq  and  Terazawa 
treatment  of  the  static  load  applied  to  the  surface  of  a  linear  elastic 
half-space.  ’  For  real  materials,  the  depth  of  influence  depends  on 
the  size  of  the  loaded  area,  the  magnitude  of  the  static  load,  and  the 
intrinsic  nonlinear  properties  of  the  material.  For  a  dynamic  load, 
there  is  also  a  finite  range  of  influence  in  real  media  where  the  dy¬ 
namic  amplitude  of  the  three-dimensional  elastic  waves  emanating  from  a 
source  goes  to  zero.  For  the  dynamic  source,  the  finite  range  of  in¬ 
fluence  is  generally  thousands  of  feet  or  often  miles;  but  for  the 
static  source,  the  finite  influence  distance  is  generally  a  few  feet  or 
inches.  Both  homogeneous  and  nonhomogeneous  materials  exhibit  a  finite 
range  of  influence,  and  in  this  report  it  is  assumed  that  the  finite 
static  depth  of  influence  combined  with  the  variation  of  the  elastic 
properties  of  a  pavement  with  depth  determines  the  magnitude  of  the  non¬ 
linear  response  of  a  pavement  to  static  and  dynamic  loads  applied  to  the 
surface. 

Only  a  small  departure  from  the  linear  stress-versus-strain  con¬ 
stitutive  equation  is  required  to  produce  a  finite  range  of  influence  of 
several  feet.  A  small  specimen  of  this  solid  would  exhibit  essentially 
linear  characteristics  during  laboratory  load-deflection  tests  or  stand¬ 
ing  wave  tests.  Only  over  a  distance  of  several  feet  in  an  actual  soil 
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formation  or  pavement  would  the  nonlinear  constitutive  equation  manifest 
itself  by  the  finite  range  of  the  stress  and  strain  field.  The  spring 
constant  for  a  real  elastic  half-space  depends  on  the  static  range  of 
influence  and  on  the  values  of  v  and  G  which  are  contained  within 
the  static  range  of  influence.  In  this  chapter,  the  foregoing  arguments 
are  placed  on  a  quantitative  basis  by  a  calculation  of  the  nonlinear 
response  of  a  one-dimensional  representation  of  an  elastic  half-space 
for  auy  distribution  of  v  and  G  with  depth. 

The  four  calculations  required  to  determine  the  three  parameters 
,  b  ,  and  e  are: 

a.  Work  done  by  the  static  elastic  deflection  of  the  pavement 
surface  beneath  the  vibrator  baseplate  (Section  3.6.3). 

b.  Static  strain  energy  in  the  pavement  (Section  3.6. 4). 

c_.  Calculation  of  average  elastic  parameters  for  a  one¬ 
dimensional  layered  system  (Sections  3.6.5  and  3.6.6). 

d.  Representation  of  the  finite  depth  of  influence  by  an  in¬ 
finite  series  (Section  3.6.7). 

3.6.3  WORK  ASSOCIATED  WITH  STATIC 
ELASTIC  DEFLECTION 


The  dynamic  load  applied  at  the  surface  is  determined  by  the 
three  parameters  k^  ,  b  ,  and  e  .  A  general  method  of  calculating 
these  three  parameters  in  terms  of  the  variation  of  v  and  G  with 
depth  is  required.  The  static  elastic  indentation  as  a  function  of 
static  load  can  then  be  calculated  from  Equation  3.1+  using  these  three 
parameters.  The  method  used  to  calculate  k^  ,  b  ,  and  e  is  an  en¬ 
ergy  method  which  equates  the  work  done  during  the  indentation  of  a 
rigid  punch  into  the  surface  of  an  elastic  half-space  to  the  increase  of 
the  potential  strain  energy  of  the  half-space.  The  work  done  by  the 
static  load  is  obtained  from  Equation  3.4  to  be 


^nnxp  ex^ 

w  =  -22-£  +  -ri  ♦ 

2  4  6 


=  W  +  W,  +  V.Vx^ 
2  e  4  e  6  e 


(3.50) 
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where  W2  ,  ,  and  Wg  are  the  coefficients  in  the  power  series  ex¬ 

pansion  of  W  .  Only  a  static  approach  to  the  calculation  of  k^  , 
b  ,  and  e  is  undertaken  in  this  report,  and  it  should  be  realized  that 
a  fully  dynamic  calculation  of  these  parameters  may  introduce  a  fre¬ 
quency  dependence. 

3.6.4  STRAIN  ENERGY  OF  AN  ELASTIC 
HALF-SPACE  WITH  A  FINITE 
DEPTH  OF  INFLUENCE 

The  elastic  strain  energy  of  a  pavement-soil  system  can  be  ex¬ 
pressed  in  terms  of  the  equilibrium  static  elastic  displacement  xg  . 

A  complete  three-dimensional  calculation  of  the  strain  energy  of  a  non- 
homogeneous  nonlinear  half-space  is  very  complicated.  The  calculation 
of  the  static  strain  energy  that  is  done  in  this  chapter  is  a  one¬ 
dimensional  approximation  of  the  actual  physical  situation.  In  this 
calculation  it  is  assumed  that  the  strain  in  the  half-space  is  in  the 
vertical  direction  only  and  is  confined  to  a  frustum  of  a  cone  whose 
upper  area  is  equal  to  the  area  of  the  vibrator  baseplate.  The  height 
of  the  frustum  of  a  cone  is  taken  to  be  equal  to  the  finite  depth  of  in¬ 
fluence  in  the  pavement,  and  the  area  of  the  lower  base  is  chosen  to 
make  the  theoretical  values  of  the  finite  depth  of  influence  agree  with 
the  experimental  values. 

The  frustum  of  the  cone  in  which  the  strain  is  assummed  to  be 
confined  is  shown  in  Figure  5-  The  volume  of  the  frustum  is  given  by 

V  =  (3.51) 

where 

V  =  volume  of  the  frustum  of  a  cone  having  a  depth  equal  to  the 
finite  depth  of  influence  in  the  pavement 

l  =  finite  depth  of  influence  of  the  static  strain  field 

4*  =  the  volume  factor  for  the  frustum  of  the  cone  which  is  cal¬ 
culated  as  follows: 


4 '  = 


1  + 


K  +  K 


2 


3 


(3.52) 
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b3 


where  «  equals  the  ratio  of  the  radius  of  the  lower  base  of  the 

frustum  to  the  radius  of  the  upper  area  of  the  frustum.  Experimental 

tests  show  that  S'  *  1+  and  k  %  3  for  AC  pavements. 

It  is  assumed  that  the  strain  energy  can  be  written  in  a  form 

analogous  to  that  for  a  perfectly  linear  elastic  system  except  that  the 

depth  of  influence  is  taken  to  be  finite.  The  elastic  strain  energy 

ll 

density  for  a  homogeneous  solid  is 

6  =  Ge2  +  (3.53) 


where 

6  =  energy  density 
e  =  strain  in  vertical  direction 
X  -  Lame  elastic  constant 
9  =  elastic  volume  dilation 

The  validity  of  Equation  3.53  for  a  material  with  a  finite  depth  of  in¬ 
fluence  arises  from  the  essentially  linear  elastic  behavior  of  most 
solids  including  PCC  and  AC  under  a  small,  uniformly  disturbed  stress. 
This  implies  that  G  and  X  have  a  physical  meaning  even  for  real 
solids.  It  is  only  over  a  long  distance  from  a  static  surface  source  on 
a  half-space  (which  produces  a  nonuniform  distribution  of  stress  and 
strain)  that  the  nonlinear  effects  manifest  themselves  by  producing  a 
definite  cutoff  range  for  the  stress  and  strain. 

An  approximation  to  the  strain  energy  density  of  a  real  half-space 
is  given  by  Equation  3.52  if  the  finite  depth  of  influence  is  introduced 
through  e  and  8  .  For  a  one-dimensional  representation  of  an  elastic 
half-space,  the  strain  and  dilation  can  be  approximated  by 
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where 


A£  =  increment  of  the  finite  depth  of  influence 

AV  »  increment  of  the  volume  of  the  frustum  of  the  cone  contain¬ 
ing  the  strain  field 

Therefore  within  the  limits  of  the  vertical  strain  approximation,  e 
and  0  are  equal,  and  the  strain  energy  density  for  a  homogeneous  half¬ 
space  can  be  written  as 


■  *  •  m 


(3.56) 


The  total  strain  energy  for  the  homogeneous  material  in  the  cylinder 


of  influence  can  be  written  as 


U  =  6V  =  ira‘ 


(3.57) 


where  U  is  the  elastic  strain  energy.  In  engineering  practice  it 
is  customary  to  use  Poisson's  ratio  in  place  of  the  Lame  elastic  con¬ 
stant.  In  terms  of  Poisson's  ratio,  Equation  3.57  Is  written  as 


2  2 
rra  QGx  ¥ 
e 


(3.58) 


where  Q  ,  the  function  of  Poisson's  ratio,  is  expressed  as 


1  -  v 
1  -  2v 


(3.59) 


For  a  nonhomogeneous  half-space  for  which  v  and  G  vary  with  depth, 


the  stx-ain  energy  is  written  as 


ira  QGx  ¥ 
e 


(3.60) 


where  QG  equals  the  value  of  QG  averaged  over  the  depth  of 
influence  l  . 


3.6.5  LAYERED  ELASTIC  MODEL  OF 
PAVEMENT  AND  SOIL  SYSTEM 


A  pavement-soil  system  generally  consists  of  a  layer  of  PCC  or  AC 


overlying  several  layers  of  crushed  rock  base  and  subbase,  all  of  which 
rest  on  the  subgrade  which  is  the  natural  soil  foundation.  A  model 
representing  the  pavement-soil  system  consists  of  a  series  of  elastic 
layers,  with  v  and  G  constant  in  each  layer  to  represent  the  pave¬ 
ment,  base,  and  subbase  overlying  an  elastic  half-space.  This  half¬ 
space  represents  the  subgrade  and  therefore  has  v  and  G  increasing 
continuously  with  depth  because  the  overburden  pressure  increases  con¬ 
tinuously  with  depth.  This  variation  of  v  and  G  with  depth  is  rep¬ 
resented  in  Figure  6. 

3.6.6  AVERAGED  ELASTIC  PARAMETERS 

The  average  value,  QG  ,  required  in  Equation  3.60  is  calculated 
by  averaging  over  the  depth  of  influence  £  as  follows: 

QG  =  Q1G1  ,  £  <  hx  (3.61) 

QG  =  j[Q1G1h1  +  Q2G2U  -  h1)J  ,  hx  <  £  <  H2  (3.62) 

*  -  t[Wl  +  Wz  *  S°3(‘  -  H2>]  •  H2  '  1  ‘  H3  (3'63) 

where  Q^G^  ,  Q2G2  ,  and  Q^G^  equal  respectively  the  value  of  QG  in 

layers  1,  2,  and  3,  h^  ,  hg  ,  and  h^  equal  respectively  the  thick¬ 
nesses  of  layers  1,  2,  and  3,  ,  Hg  ,  and  equal  the  sums  of 

pavement  layer  thicknesses,  and  where 

Hg  =  h1  +  h2  (3.61) 

H3  =  h1  +  h2  +  h3  (3.65) 

Equation  3.6l  is  valid  if  the  depth  of  influence  extends  only  into  the 
upper  layer,  while  Equation  3.62  is  valid  if  £  extends  into  the  second 

layer,  and  finally  Equation  3.63  is  valid  if  £  extends  into  the  third 

layer.  The  generalization  of  Equations  3.61-3.63  to  more  than  three 

layers  is  obvious.  Equations  3.61-3.63  are  valid  only  if  the  depth  of 


J'6 
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Figure  6.  Layered  model  of  pavement 


influence  is  contained  within  the  pavement,  base,  and  subbase,  i.e.,  the 
layers  in  which  v  and  G  are  constant.  If  the  depth  of  influence  ex¬ 
tends  into  the  subgrade,  then  the  proper  average  is 


QG 


T  [VA  ♦  Wz  +  %°3h3  +  Qs°= 


£  >  h3  (3.66) 


where  G  equals  the  shear  modulus  of  the  subgrade,  Q  equals  the 

S  S 

value  of  Q  given  by  Equation  3.59  but  for  the  subgrade,  and  Q  G 

s  s 

equals  the  average  value  of  Q  G  in  the  interval  £  -  H_  .  The 

S  S  j 

product  Q  G  for  soil  at  a  depth  z  below  the  surface  of  the  pavement 
s  s 

can  be  written  as  a  Taylor  series  expansion  as  follows: 

(z  -  Hj 


0  0  d(QsGs) 

Q  (z)G  (z)  =  Q°G°  +  - 

s  s  s  s  dz 
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2  dz2 


(z  -  h 3r  +. 


(3.67) 


where  QGG°  equals  the  product  Q  G  for  the  soil  at  the  very  top  of 
s  s  s  s 

the  subgrade,  i.e.,  the  value  of  Q  G  at  £  =  H_  .  The  average  value 

S  S  j 

Q  G  appearing  in  Equation  3.66  is  assumed  to  be  given  as 
s  s  % 


Q  G 
s  s 


=  1  -1h3  / 


}  G  dz 
s  s 


(3.68) 


Using  Equation  3.67,  the  average  becomes 


rs-  =  q°g°  *  i 
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and  this  value  is  used  in  Equation  3.66. 


1*8 


3.6.7  SERIES  EXPANSION  OF  THE  FINITE 
DEPTH  OF  INFLUENCE 

The  physics  of  pavements  and  soils  enters  through  a  calculation 
of  the  depth  of  influence,  l  ,  from  a  fundamental  knowledge  of  the 
molecular  structure  of  pavement  and  soil  materials  and  from  a  knowledge 
of  the  cracks,  voids,  dislocations,  and  other  microscopic  flaws  in  these 
materials.  A  basic  calculation  of  l  from  a  theoretical  physics  point 
of  view  has  not  been  carried  out  in  this  study.  In  this  report,  the 
static  depth  of  influence  is  assumed  to  be  given  by 

l  =  +...  (3.70) 

where  ,  and  are  coefficients  of  the  power  series  expan¬ 

sion  of  the  finite  depth  of  influence,  and  depend  on  the  radius  of  the 
contact  area  of  the  load  with  the  pavement  surface  and  on  the  values  of 
v  and  G  of  the  pavement-soil  system.  The  choice  of  the  form  of  Equa¬ 
tion  3.70  is  dictated  by  the  even  powers  of  the  work  function  in  Equa¬ 
tion  3.50.  A  rigorous  calculation  of  the  values  of  ,  and  Z^ 

has  not  been  accomplished  theoretically.  In  this  report  the  values  of 
,  Z^  ,  and  ^  are  estimated  from  the  experimental  values  of  k^  , 
b  ,  and  e  . 

3.7  CALCULATION  OF  THE  NONLINEAR 
PARAMETERS 

All  of  the  quantities  required  for  the  determination  of  k^  , 
b  ,  and  e  have  been  calculated  in  the  previous  five  sections.  Now  the 
energy  principle,  that  the  static  potential  strain  energy  must  equal  the 
work  done  during  the  static  elastic  deflection  of  the  vibrator  into  the 
pavement  surface,  is  applied  to  the  one-dimensional  pavement  model. 

This  condition  gives  the  connection  between  the  parameters  k^  ,  b  , 
and  e  and  the  corresponding  elastic  and  inelastic  pavement  parameters. 
In  this  section  the  parameters  k  t  b  ,  and  e  are  calculated  for  a 
pavement  system  consisting  of  a  series  of  layers  overlying  a  subgrade. 


3.7.1  GENERAL  EXPRESSION  FOR  THE 
PAVEMENT  PARAMETERS 

The  potential  strain  energy  is  calculated  using  Equations  3.60- 
3.70  and  can  he  put  in  the  form: 

U  =  U2xJ  ♦  V*  ♦  B(j*«  (3.71) 

where  U2  ,  ,  and  Ug  are  the  coefficients  of  the  power  series  ex¬ 

pansion  of  the  elastic  strain  energy.  The  particular  mathematical  form 
for  the  potential  strain  energy  given  by  Equation  3.71  is  chosen  in 
order  to  agree  with  the  form  of  the  work  function  in  Equation  3.50.  The 
values  of  ,  and  Ug  are  calculated  by  substituting  the  value 

of  l  given  by  Equation  3-70  into  the  expression  QG  in  Equation  3.60 

2 

and  expanding  in  a  power  series  in  xg  .  The  values  of  ,  and 

Ug  will  depend  on  the  magnitude  of  the  depth  of  influence,  which  is 
determined  by  ,  and  ,  and  on  the  v  and  G  values  within 

the  depth  of  influence.  The  coefficients  k^  ,  b  ,  and  e  are  ob¬ 


(3.72) 


(3.73) 


In  the  following  sections,  the  general  expressions  for  the  param¬ 
eters  k^  ,  b  ,  and  e  given  by  Equation  3.73  are  evaluated  for 
specific  cases  of  layered  systems  and  layered  systems  overlying  a 


tained  by  setting  W  =  U  or  equivalently 


W2  =  U2 


from  which  it  follows  that 


W6  =  U6 


k00  =  2U2 
b  =  kuh 

e  =  6U, 


The  values  of  kQ0  ,  b  ,  and  e  will  depend  on  ,  and 
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subgrade  whose  shear  modulus  increases  uniformly  with  depth. 

3.7.2  ONE  ELASTIC  LAYER 

For  £  <  h^  ,  the  stress  and  strain  extend  only  into  the  upper 
pavement  layer.  Equations  3.60,  3.6l,  and  3.70  give: 


2wa  Q^G  ¥ 


00 


(3. Ik) 


b  =  - 


(3.75) 


e  = 


6na2Q1G154' 


(3.76) 


where  6  is  the  function  of. the  expansion  coefficients  of  the  finite 
depth  of  influence  and 
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(3.77) 


The  results  of  Equations  3.74-^. 77  are  valid  for  the  homogeneous  half¬ 
space. 

3.7.3  TWO  ELASTIC  LAYERS 

For  h^  <  l  <  Hg  ,  the  stress  and  strain  extend  into  the  second 
layer  and  Equations  3.60,  3-62,  and  3.70  give: 
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where  p  is  the  function  of  the  expansion  coefficients  of  she  finite 
depth  of  influence  and 


P  = 


(3.81) 


Equations  3. 78-3.80  reduce  to  Equations  3. 7^-3. 76  when  ,  as 

they  should  because  this  condition  reproduces  the  homogeneous  half-space. 

3,-l.h  THREE- LAYER  SYSTEM  AND  THE 
GENERALIZATION  TO  n  LAYERS 

For  the  three-layer  case  with  <  l  <  H^  corresponding  to  the 
case  where  the  stress  and  strain  extend  into  the  third  layer.  Equations 
3.60,  3.63,  and  3.70  give: 
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+  Q3G3^0V  (3.83) 


Q0GJ  +  h0{ Q0G  -  Q_G,)  +  6Q0G0£n )  (3.8M 


The  generalization  to  the  case  of  more  than  three  layers  is  simple. 

Consider  a  system  of  n  layers  where  H  <  &  <  H  and  which  corre- 

n-1  n 

sponds  to  the  case  in  which  the  finite  influence  depth  extends  into  the 
nth  layer.  For  this  case  the  coefficients  are 


2 

2va  V 


00  2 
0 


n-1 

Zh.(Q.C.  -  Q  G  )  +  Q  G  l. 

111  n  n  t  n  C 


b  =  - 


i=l 
n-1 


2  Y  h.  (Q.C.  -  Q  G  )  +  Q  G  t 
^111  n  n  n  n  C 


i=l 


8rO 


(3.86) 


52 


e  = 


6ira2Y 


n-1 


p  S  h. (Q. G.  -  Q  G  )  +  6Q G  l 
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(3.87) 


0  L  i=l 
3.7.5  LAYERED  SYSTEM  WITH  SUBGRADE 

If  the  n  layered  pavement  overlies  a  subgrade  and  £  >  Hn  so 
that  the  stress  and  strain  extend  into  the  soil  medium,  then  Equations 
3.60  and  3.66-3.70  give 
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where  ,  Eg  ,  Kg  ,  ,  and  <g  are  the  expansion  coeffi¬ 

cients  appearing  in  the  expressions  for  k^,  b  ,  and  e  and 
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For  a  soil  formation  by  itself,  lu  =  0  and  Hn  =  0  so  that 
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where  QGGG  depends  on  the  values  of  v  ,  Poisson's  ratio  for  the 
s  s  s 

subgrade,  and  at  the  soil  surface. 
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3.7.6  ALGEBRAIC  SIGNS  OF  THE  NON¬ 
LINEAR  PARAMETERS  b  AND  e 

The  minus  sign  appearing  in  the  expression  for  the  nonlinear  co¬ 
efficient  b  arises  from  the  inverse  dependence  of  U  on  l  .  From 
Equation  3.75  it  follows  that  b  <  0  for  a  homogeneous  half-space.  For 
a  nonhomogeneous  half-space,  the  sign  of  the  coefficient  b  depends  on 
the  variation  of  v  and  G  with  depth  in  the  pavement  and  in  the  soil. 
Consider  first  the  case  where  the  depth  of  influence  extends  only  into 
the  pavement,  Z  <  ,  and  Equation  3.86  is  valid.  Then  Equation  3.86 

shows  that  if  the  values  of  Q.G^  decrease  with  depth,  as  is  the  gen¬ 
eral  case  with  pavements,  then  b  <  0  .  If  Q^G^  increases  with  depth, 
then  b  <  0  or  b  >  0  ,  depending  on  the  rate  of  increase  of 
with  depth.  If  the  rate  of  increase  of  Q^G^  with  depth  is  suffi¬ 
ciently  rapid,  then  there  is  a  possibility  that  b  >  0  .  The  same  re¬ 
sults  are  valid  when  the  depth  of  influence  extends  into  the  subgrade. 
Most  cases  of  interest  to  soil  and  pavement  engineers  will  have  b  <  0  . 

3.7.7  CALCULATION  OF  THE  FINITE  DEPTH 
OF  INFLUENCE  AND  THE  DETERMINA¬ 
TION  OF  LAYER  THICKNESS 

It  has  been  shown  in  Sections  3. 7- 2-3. 7- 5  that  different  expres¬ 
sions  for  the  coefficients  kQ0  ,  b  ,  and  e  must  be  chosen  according 
to  the  magnitude  of  the  depth  of  influence,  i.e.,  depending  on  the  par¬ 
ticular  pavement  layer  into  which  the  static  depth  of  influence  extends. 
The  static  finite  depth  of  influence  is  expected  to  depend  on  the  base¬ 
plate  radius  of  the  vibrator  and  the  static  load  that  the  vibrator  apr 
plies  to  the  pavement  surface.  Therefore  different  expressions  for  the 
coefficients  ,  b  ,  and  e  must  be  used  if  a  series  of  baseplate 

radii  are  selected  for  the  vibrator.  Equations  3. 7^-3. 99  show  that  the 
change  in  the  equations  describing  k^  ,  b  ,  and  e  occurs  when 

,  i.e.  ,  when  Z^  passes  through  the  successive  interfaces  of 
the  pavement  system.  The  values  of  the  baseplate  radius  for  which 

Zn  =  H.  will  be  called  the  critical  radii.  Therefore  the  critical  radii 
0  i 

are  defined  by 
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(3.100) 


Vaci>  =  Hi 


where  a  ,  ,  a  „  ,...a  .  are  the  critical  contact  radii  of  the  vibra- 
cl  *  c2  *  ci 

tor  baseplate  corresponding  to  the  interface  at  the  depths  H.  . 

The  zero^  -order  coefficient,  £q  ,  of  the  static  finite  depth  of 
influence  can  be  calculated  in  terms  of  the  radius  of  the  vibrator  base¬ 
plate,  the  critical  radii  for  the  layered  system,  and  the  elastic  con¬ 
stants  of  each  layer  of  the  system.  The  value  of  £,  for  a  homogeneous 
elastic  half-space  can  be  obtained  by  equating  the  value  of  given 

by  Equation  3-71*  with  the  value  of  the  spring  constant  for  a  linear 
elastic  half-space,  with  the  result: 

lQ  =  |  naYB  (3.101) 


where  B  is  a  function  of  Poisson's  ratio  for  the  homogeneous  half¬ 
space  and  is  defined  by: 
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Equation  3.101  shows  that  is  proportional  to  the  radius  of  the 

contact  area  of  the  vibrator  which  is  exerting  the  static  load  on  the 
surface  of  the  half-space. 

For  a  layered  pavement  system,  is  given  by  the  following 

expressions: 
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where  the  values  of  B  for  each  pavement  layer  are  given  by: 


(3.107) 


and  v.^  ,  are  the  values  of  Poisson's  ratio  for  the  1,  2,...i 

layers  of  the  pavement. 

In  Section  ^.1  it  will  be  shown  that  the  critical  radii  of  a 


pavement  system  can  be  measured  experimentally  using  a  mechanical  vi¬ 
brator,  so  that  in  principle  the  function  £g(a,ac^)  can  be  evaluated 
The  value  of  the  first  critical  radius  is  related  to  the  thickness  of 


the  first  layer  in  the  following  manner: 
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In  a  similar  manner  it  is  easy  to  show  from  Equations  3.103-3.105  that 
the  thicknesses  of  the  second  and  third  layers  are  given  by: 


-  ?  •'Vv  -  Vl> 

(3.109) 

(3.110) 

Therefore  the  thickness  of  each  pavement  layer  can  be  determined  if  the 
values  of  the  critical  radii  can  be  measured. 

The  values  of  £„  and  £,  can  be  obtained  in  terms  of  the 
values  of  ,  b  ,  and  e  using  Equations  3.75,  3.76,  3.79,  3.80, 
3.83,  3.81+,  3.86,  3.87,  3.89,  and  3.90.  Therefore  measurements  of  b  , 
e  ,  and  the  critical  radii  are  necessary  for  the  complete  determination 
of  the  finite  depth  of  influence.  A  theoretical  calculation  of  lQ  , 

J ,  and  would  probably  involve  either  a  molecular  theory  to  de¬ 

scribe  the  nonlinear  behavior  of  the  molecular  bonds  in  AC  and  PCC  or  a 
semi  empirical  nonlinear  elastic  theory,  which  has  not  been  considered  in 
this  report.  Numerical  values  of  ZQ  ,  ,  and  ^  can  be  obtained 

from  measured  values  of  k^  ,  b  ,  and  e  ,  provided  that  the  elastic 
constants  and  the  layer  thicknesses  of  the  pavement  are  known.  The 
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numerical  values  of  the  coefficients  l. 


J&2  ,  and 


that  are  de¬ 


scribed  in  Section  4.2.2  were  estimated  in  this  manner. 

3.7-8  EFFECT  OF  BASEPLATE  RADIUS 

b  ,  and  e  depend  on  the  following 


k00  ’ 


The  coefficients 
quantities : 

a.  The  elastic  parameters  and  of  the  pavement  layers. 


b.  The  elastic  parameters 

c.  The  coefficients  l 


v  and  G  of  the  subgrade . 
s  s 


The  coefficients  £, 
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loaded  area  on  the  surface,  and  therefore  k 

The  values  of  k 
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departure  from  the  linear  stress-strain  constitutive  equation. 


pend  on  this  radius. 


00 


00  *  ’ 
b  ,  and  e 


and  e  also  de¬ 
can  be  obtained 


as  a  function  of  baseplate  radius  by  combining  Equations  3.74-3.87  with 
the  equations  giving  2^  as  a  function  of  baseplate  radius.  The  basic 
forms  of  the  parameters  ,  b  ,  and  e  and  the  dynamic  stiffness, 

S  ,  are  given  in  Figures  7»  8,  9,  and  10,  respectively,  for  the  cases 
where  QG  decreases  with  depth  (as  in  pavements)  and  for  the  cases 
where  QG  increases  with  depth.  The  measured  degree  of  nonlinearity  of 
a  pavement-soil  system  should  depend  on  the  radius  of  the  vibrator  base¬ 
plate.  For  very  large  values  of  the  radius,  the  slopes  of  the  param¬ 
eters  k^  ,  b  ,  and  e  with  respect  to  the  baseplate  radius  are  de¬ 
termined  by  the  deeper  pavement  or  soil  layers,  whereas  the  values  of 
these  three  parameters  are  determined  by  the  elastic  properties  of  all 
of  the  layers  which  are  contained  in  the  depth  of  influence  2^  .  For 
very  small  values  of  the  baseplate  radius,  the  slopes  of  kQQ  ,  b  ,  and 
e  with  respect  to  the  baseplate  radius  are  determined  by  the  elastic 
properties  of  the  surface  and  near-surface  layers.  The  parameters  b 
and  e  are  discontinuous  functions  of  the  radius  of  the  loaded  area; 
the  discontinuities  occur  when 

two  layers,  i.e.,  when  2,  =  H.  where  H.  is  the  total  thicKness  or  1 

layers.  The  value  of  the  discontinuities  at  x,„  =  h.  are 


passes  through  the  boundary  between 

H.  is  the  total  thickness  of 
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b.  CASE  IN  WHICH  SHEAR  MODULUS  INCREASES  WITH  DEPTH  (SOILS) 


Figure  7.  Theoretical  dependence  of  the  spring  constant,  k  , 
on  the  radius  of  the  vibrator  baseplate  ^ 
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b.  CASE  IN  WHICH  SHEAR  MODULUS  INCREASES  WITH  DEPTH  (SOILS) 


Figure  8.  Theoretical  dependence  of  the  third-order  nonlinear 
parameter,  b  ,  on  the  radius  of  the  vibrator  baseplate 
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b.  CASE  IN  WHICH  SHEAR  MODULUS  INCREASES  WITH  DEPTH  (SOILS) 


Figure  9- 
parameter, 


Theoretical  dependence  of  the  fifth-order  nonlinear 
e  ,  on  the  radius  of  the  vibrator  baseplate 
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b.  CASE  IN  WHICH  SHEAR  MODULUS  INCREASES  WITH  DEPTH  (SOILS) 

Figure  10.  Theoretical  dependence  of  the  dynamic  stiffness, 
on  the  radius  of  the  vibrator  baseplate 
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ae  =  (P  -  6)  H.(4i01  -  Qln0ln)  (3.112) 

*o 

The  values  of  the  discontinuities  in  b  and  e  can  be  positive  or 
negative  depending  on  whether  QG  is  a  decreasing  or  an  increasing 
function  of  depth,  as  indicated  in  Figures  8  and  9- 

3.7-9  DETERMINATION  OF  THE  SHEAR  MOD¬ 
ULI  OF  THE  PAVEMENT  LAYERS 

The  shear  modulus  of  each  pavement  layer  can  be  obtained  from  the 
values  of  k^  for  a  series  of  baseplate  radii.  Computer  programs  have 
been  developed  for  calculating  the  values  of  from  the  measured 

dynamic  load-deflection  curves.  Therefore,  if  dynamic  load-deflection 
curves  are  measured  for  a  series  of  vibrator  baseplate  radii,  it  is 
possible  to  obtain  kQ0  as  a  function  of  baseplate  radius.  The  elastic 
shear  modulus  of  each  pavement  layer  is  expressed  in  terms  of  V  by 
using  Equations  3. 7^-3. 89  with  the  result: 
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These  expressions  for  the  shear  moduli  can  be  expressed  in  terms  of  the 
baseplate  radius  and  the  critical  baseplate  radii  by  using  Equations 
3.103-3.110,  with  the  result: 
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where  0o(a)  ,  0^(a)  ,  and  0^(a)  are  functions  of  the  baseplate 

radius  and  Poisson's  ratio  of  the  successive  pavement  layers.  These 
functions  are  given  by 
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Therefore  the  values  of  the  shear  modulus  of  each  layer  of  the 
pavement  system  can  be  calculated  if  the  values  of  Poisson's  ratio  are 
assumed  to  be  known  for  each  pavement  layer  and  if  the  function  k^^a) 
is  known  for  a  series  of  values  of  the  baseplate  radius.  This  series 
must  include  a  baseplate  radius  sufficiently  large  to  make  the  static 
finite  depth  of  influence  penetrate  into  the  deepest  layer  whose  shear 
modulus  is  to  be  determined.  The  function  k^a)  can  be  obtained 
from  measured  values  of  the  dynamic  stiffness  S(a,Fs,Fd)  .  The  shear 
modulus  of  the  subgrade  can  be  determined  by  this  method  if  a  suffi¬ 
ciently  large  vibrator  baseplate  is  used.  It  is  clear  also  from  Equa¬ 
tions  3.97-3-99  that  this  method  is  also  applicable  to  soil  formations 
and  may  prove  to  be  a  useful  method  of  shallow  f  -.bsurface  geophysical 
exploration. 

The  unique  feature  of  this  method  of  determining  the  variation 
of  the  shear  modulus  with  depth  is  the  one-to-one  correspondence  between 
the  vibrator  baseplate  radius  and  the  finite  depth  of  influence  in  the 
pavement.  Each  pavement  layer  can  be  detected  individually  by  choosing 
a  proper  baseplate  size.  This  one-to-one  correspondence  does  not  exist 
in  the  wave  propagation  method  of  determining  subsurface  structure,  be¬ 
cause  in  this  case  a  single  point  on  a  measured  Rayleigh  wave  dispersion 
curve  is  influenced  by  the  values  of  the  shear  modulus  at  all  depths  be¬ 
neath  the  surface. 

The  basic  purpose  of  applying  pavement  layers  over  a  subgrade  is 
to  distribute  the  load  which  is  applied  to  the  pavement  surface  in  such 
a  way  as  to  protect  the  comparatively  weak  subgrade.  A  measure  of  the 
load  distribution  characteristics  of  a  pavement  is  given  by  the  angle  of 
the  frustum  of  the  cone  in  which  the  strain  (and  stress)  in  the  pave¬ 
ment  system  is  assumed  to  be  confined.  This  angle  is  shown  in  Figure  5 
and  is  defined  as  follows: 


tan  0 


-  1 

D  =  (k  -  1  )a 


(<  -  ill 


(3.1214) 


where  0^  is  the  angle  of  stress  distribution.  For  a  homogeneous  elas¬ 
tic  half-space  this  angle  is  given  by 
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(3.125) 
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Similar  expressions  for  the  layered  pav  jment  system  can  be  obtained  by 
using  Equations  3.104-3.106.  For  the  homogeneous  elastic  half-space, 

0p  is  independent  of  the  baseplate  radius;  but  for  layered  systems, 
this  angle  does  depend  somewhat  on  the  baseplate  radius.  The  values  of 
k  are  determined  experimentally  by  requiring  that  Equations  3.108-3.110 
correctly  predict  the  pavement  layer  thicknesses.  This  requirement 
gives  a  value  of  ¥  from  which  the  value  of  k  is  immediately  obtained 
by  using  Equation  3.52.  A  typical  value  of  for  pavements  is 

0D  »  70  deg. 

3.8  EFFECTS  OF  THE  MECHANICAL  CHARAC¬ 
TERISTICS  OF  THE  VIBRATOR  ON  THE 
MEASURED  VALUES  OF  DYNAMIC  STIFFNESS 

3.8.1  VIBRATOR  AND  PAVEMENT  PARAMETERS 
INFLUENCING  THE  DYNAMIC  STIFFNESS 

The  dynamic  stiffness  of  a  pavement  is  measured  by  using  a  me¬ 
chanical  vibrator  to  determine  the  load-deflection  curves.  The  shape 
of  these  curves  depends  on  the  values  of  the  nonlinear  parameters  b 
and  e  .  The  experimental  results  of  Section  4.1  and  the  theoretical 
results  of  Section  3.5  indicate  that  the  nonlinear  parameters  b  and  e 
are  not  zero  for  pavements.  In  this  case.  Equations  3.22-3.28  show  that 
in  addition  to  the  quantities  listed  in  Section  2.4  pertaining  to  the 
linear  oscillator,  the  dynamic  stiffness  of  a  pavement  measured  by  a 
vibrator  also  depends  on  the  following  quantities: 

a.  Two  nonlinear  elastic  parameters  describing  the  pavement. 

b.  The  static  load  (or  static  equilibrium  displacement). 

c_.  The  dynamic  force  generated  by  the  vibrator. 

d..  The  baseplate  radius  of  the  vibrator. 

Because  various  vibrators  have  different  values  of  F  and  F  ,  it 

3  l) 

follows  that  the  dynamic  stiffness  value  measured  by  different  vibrators 
at  the  same  locality  will  in  general  be  different.  The  dynamic  stiff¬ 
ness  measured  by  a  vibrator  thus  depends  on  the  physical  cnaracterist' cs 
of  the  vibrator  used  for  the  measurements.  For  flexible  pavement  it  is 
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in  part  the  dependence  of  the  dynamic  stiffness  on  that  produces 
the  discrepancy  in  the  measured  values  of  dynamic  stiffness  using  two 
different  vibrators  with  different  values  of  at  the  same  pavement 

location.  All  other  factors  being  equal,  two  vibrators  operating  at 
different  levels  of  dynamic  force  will  not  yield  the  same  value  of  dy¬ 
namic  stiffness.  The  lower  value  will  occur  in  the  region  of  large 
dynamic  load. 

3.8.2  THEORY  OF  THE  ANOMALOUS  BEHAVIOR 
OF  THE  DYNAMIC  STIFFNESS 

The  spring  constant  expression  shows  that  the  linear  theory  of 
elasticity  implies  a  value  of  the  dynamic  stiffness  which  is  essentially 
proportional  to  the  baseplate  radius  of  the  vibrator  mass.  Experimental 
results  presented  in  Section  h.l  contradict  this  basic  conclusion  of 
linear  elasticity.  Some  of  the  experimental  data  in  Section  h.l  indi¬ 
cate  that  for  the  same  vibrator  (same  static  and  dynamic  loads),  the 
measured  values  of  dynamic  stiffness  are  roughly  independent  of  base¬ 
plate  radius.  In  fact,  some  of  the  data  in  Section  l*.l  show  that  for 
different  vibrators  (different  static  loads),  it  is  sometimes  possible 
to  measure  smaller  values  of  the  dynamic  stiffness  with  a  large  base¬ 
plate  than  with  a  small  baseplate  at  the  same  location.  This  conflict 
between  the  experimental  results  and  the  predictions  of  linear  elastic 
theory  will  be  referred  to  as  the  anomalous  behavior  of  the  dynamic 
stiffness. 

The  anomalous  behavior  of  the  dynamic  stiffness  can  be  explained 
in  terms  of  the  nonlinear  response  of  a  pavement  to  a  static  and  dynamic 
load.  The  anomalous  result,  that  the  dynamic  stiffness  does  not  appear 
to  var;  directly  with  the  radius  of  the  baseplate  (as  it  does  for  a 
linear  homogeneous  half-space),  is  due  to  the  finite  depth  of  influence 
which  makes  only  a  very  slowly  increasing  function  of  baseplate 

radius  except  for  really  small  radii  (see  Figure  7a).  This  is  espe¬ 
cially  true  for  the  case  of  a  concrete  pavement  where  a  very  strong 
upper  layer  overlies  considerably  weaker  layers  of  base  and  subbase.  In 
this  case,  the  value  of  k^  is  determined  primarily  by  the  strong 
upper  layer,  and  the  major  part  of  this  value  can  be  obtained  with  a 
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very  small  baseplate  radius  as  shown  in  Figure  la.  The  weaker  lower 
layers  tend  to  flatten  out  the  k^-versus-baseplate-radius  curve. 

Thus,  measurements  of  k^  using  different  sized  baseplates  would  yield 
values  of  k^  considerably  less  different  than  they  would  be  if  the 
pavement  were  a  homogeneous  half-space  with  elastic  parameters  equal  to 
those  of  the  upper  layer  of  the  layered  system.  According  to  the  non¬ 
linear  theory,  k^  is  a  very  slowly  increasing  function  of  baseplate 
radius  for  baseplate  radii  larger  than  the  first  critical  radius.  There¬ 
fore,  for  sufficiently  large  baseplate  radii,  the  values  of  dynamic 
stiffness  of  rigid  pavements  measured  for  fixed  static  and  dynamic  loads 
should  be  essentially  independent  of  the  baseplate  radius. 

The  fact  that  the  values  of  k^^  are  similar  for  different  base¬ 
plate  radii  which  are  sufficiently  large  implies  that  the  values  of  two 

k  (F  ,a)  functions  corresponding  to  two  different  baseplate  sizes  are 
u  s 

also  similar  for  some  pavements.  In  this  case  it  is  possible,  due  to  a 

local  minimum  of  the  function  krt(F  ,a)  for  some  value  of  F  ,  for 

kQ(Fs,a)  for  a  small  baseplate  radius  and  small  static  load  to  be 

larger  than  kQ(Fs,a)  for  a  large  radius  and  a  large  static  load  as 

shown  in  the  rigid  pavement  curve  in  Figure  ha.  If  the  values  of  k^Q 

are  very  different  for  two  baseplate  radii,  the  corresponding  two  values 

of  the  function  kQ(Fs,a)  will  also  be  very  different.  In  particular, 

the  function  k  (F  ,a)  for  the  larger  radius  will  be  larger  than  the 
u  s 

corresponding  function  for  the  smaller  radius  for  all  values  of  F  ,  as 

s 

shown  in  the  flexible  pavement  curves  of  Figure  ha.  This  case  would 
then  represent  the  normal  situation  in  which  large  baseplates  produce 
slightly  larger  values  of  dynamic  stiffness  than  do  smaller  baseplates. 
This  normal  situation  is  expected  to  occur  for  weak  pavements  where 
there  is  not  a  great  difference  between  the  elastic  parameters  of  the 
upper  and  lower  pavement  layers.  The  anomalous  behavior  of  the  dynamic 
stiffness  measurements  on  pavements  can  be  explained  in  terms  of  the 
relative  constancy  of  the  function  k^^a)  ,  which  implies  that  the 
values  of  the  dynamic  stiffness  measured  on  rigid  pavements  by  vibrators 
with  large  and  small  baseplates  will  be  roughly  the  same.  On  the  other 
hand,  linear  elastic  theory  implies  that  the  value  of  the  dynamic 
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stiffness  should  be  directly  proportional  to  the  baseplate  radius.  The 
experimental  results,  then,  contradict  some  of  the  predictions  of  the 
linear  elastic  theory. 
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4.  EVALUATION  OF  THEORETICAL  RESULTS  OF  THE 
DYNAMIC  STIFFNESS  STUDY 

l+.l  EXPERIMENTAL  PROGRAM 

The  basic  goal  of  the  experimental  program  was  to  provide  data 
to  establish  the  use  of  the  nondestructive  evaluation  method  for  airport 
pavements.  Nondestructive  pavement  evaluation  methodology  based  on  cor¬ 
relations  between  nondestructive  data  and  existing  Federal  Aviation  Ad¬ 
ministration  evaluation  criteria  is  given  in  Volume  I  of  this  report. ^ 
The  specific  objectives  of  the  experimental  program  were: 

a.  Determination  of  the  dynamic  stiffness  as  a  function  of  fre¬ 
quency,  dynamic  load,  static  load,  and  vibrator  baseplate 
radius . 

b.  Validation  of  the  theoretical  procedures  and  predictions  of 
the  nonlinear  mechanical  model  of  pavements. 

If.  1.1  MECHANICAL  VIBRATORS  USED  FOR 
NONDESTRUCTIVE  TESTING  OF 
PAVEMENTS 

Dynamic  stiffness  data  were  obtained  from  several  vibrators 
including : 

a.  The  WES  l6-kip  hydraulic  vibrator 

b.  The  Dynaflect 

£.  The  Road  Rater 

d_.  The  Civil  Engineering  Research  Facility  (CERF)  vibrator 

e_.  The  WES  9-kip  eccentric  mass  vibrator 
Each  vibrator  has  its  own  basic  dynamic  characteristics  including  weight, 
maximum  dynamic  force,  and  baseplate  radius.  A  comparison  of  these  four 
vibrators  is  made  in  Table  1.  According  to  the  nonlinear  vibration 
theory  of  pavements  developed  in  Section  3.5,  the  dynamic  characteris¬ 
tics  of  a  vibrator  will  be  reflected  in  the  value  of  the  dynamic 
stiffness  measured  on  a  given  pavement. 

The  experimental  measurements  of  the  dynamic  load-deflection 
curves  measured  by  WES  were  made  primarily  at  a  frequency  of  15  Hz,  but 
also  in  the  range  of  10-25  Hz.  Most  of  the  data  taken  were  for  a  static 
load  of  1 6  kips  and  for  a  dynamic  load  which  varied  from  5-ll+  kips. 

Some  dynamic  stiffness  measurements  were  performed  using  a  different 
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Table  1 


Mechanical  Characteristics  of  Vibrators 


Static 

Maximum 

Load 

Dynamic  Load 

Vibrator 

kips 

kips 

WES  l6-kip 

16 

15 

Dynaflect 

1.6 

0.5 

Road  Rater 

400 

1.5 

0.75 

505 

1.5 

0.75 

510 

2.0 

1.5 

550 

4.0 

3.0 

CERF 

6.75 

5 

WES  9-kip 

9 

8 

Effective 
Baseplate 
Radius,  in. 

Contact 
Area 
pq  in. 

Operating 

Frequency 

Hz 

9 

245 

5-100 

*0.9  (concrete)* 

2.54 

8 

*3.2  (asphaltic 
concrete 

32.0 

20-30 

4.22** 

56 

20-30 

10-40 

10-40 

6 

113.1 

15,25 

9.5 

283.5 

5-60 

*  Road  Rater  applies  load  through  two  4-  by  7-in.  rectangular  steel 
pads  spaced  6  in.  apart. 

**  Dynaflect  applies  load  through  two  4-in.  wide  and  l6-in.-diam  steel 
wheels. 
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vibrator  arrangement  which  allowed  for  a  variation  of  the  static  load 
in  the  range  of  5-50  kips  by  the  addition  of  a  series  of  lead  weights 
to  the  vibrator  load. 

4.1.2  EFFECTS  OF  DYNAMIC  LOAD 

Typical  dynamic  load-deflection  curves  for  AC  pavements  appear 
in  Figure  11,  while  Figure  12  gives  these  curves  for  PCC  pavements. 

These  load-deflection  curves  correspond  to  a  frequency  of  15  Hz  and  a 
static  load  of  1 6  kips.  The  dynamic  stiffness  is  the  ratio  of  the  dy¬ 
namic  load  to  the  corresponding  value  of  the  dynamic  deflection,  and  can 
be  evaluated  for  each  point  on  the  dynamic  load-deflection  curve.  The 
value  of  the  dynamic  stiffness  corresponding  to  F^  =  0  has  been  desig¬ 
nated  by  Sq  in  Equation  3.21  and  can  be  obtained  from  Figures  11  and 
12  by  measuring  the  slope  of  the  load-deflection  curves  at  the  origin. 

The  dynamic  stiffness  values  corresponding  to  the  dynamic  load- 
deflection  curves  presented  in  Figures  11  and  12  are  given  in  Figures 
13a  and  b.  For  the  WES  l6-kip  vibrator  operating  on  AC  pavements,  the 
dynamic  stiffness  is  generally  a  decreasing  function  of  the  dynamic  load, 
while  for  rigid  pavements  the  dynamic  stiffness  is  essentially  indepen¬ 
dent  of  the  dynamic  load.  However,  the  10-Hz  dynamic  load-deflection 
curves  shown  in  Figure  l4  indicate  that  it  is  possible  under  certain 
conditions  for  the  dynamic  stiffness  to  be  an  increasing  function  of  the 
dynamic  load  and  to  have  a  local  minimum  value  at  a  specific  value  of 
the  dynamic  load.  The  results  in  Figures  11  and  12  show  that  stiff 
pavements  have  dynamic  load-deflection  curves  which  are  more  linear  than 
the  dynamic  load-deflection  curves  for  the  more  flexible  pavements. 

4.1.3  CRITICAL  FREQUENCIES 

The  dynamic  load-deflection  curves  for  a  given  pavement  have  been 
measured  for  a  series  of  frequencies,  from  which  an  estimate  of  the 
critical  frequency  for  a  vibrator-pavement  system  can  be  made.  The 
critical  frequency  has  been  defined  as  the  frequency  for  which  the  dy¬ 
namic  load-deflection  curves  are  nearly  linear  for  small  values  of  . 
Figures  l4  and  15  indicate  that  the  critical  frequency  of  the  WES  l6-kip 
vibrator  operating  on  medium  strength  AC  pavements  is  about  15.5  Hz, 
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Figure  12.  Dynamic  load-deflection  curves  for  different  locations  on  PCC  pavements 
showing  extreme  linearity  of  curves 


0  2  4  6  8  10  12  14 

DYNAMIC  LOAD,  KIPS 
a.  FLEXIBLE  PAVEMENTS 


0  2  4  6  8  10  12  14 


DYNAMIC  LOAD,  KIPS 
b.  RIGID  PAVEMENTS 

Figure  13.  Measured  values  of  dynamic  stiffness  versus 
dynamic  load  ( 15  Hz) 
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Figure  lit.  Possible  experimental  evidence  for  the  existence  of  a 
critical  frequency  value  of  *15  Hz  for  AC  pavement 
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while  other  data  indicate  that  for  PCC  pavements  it  is  about  lL.5  Hz. 
Estimates  of  the  critical  frequency,  resonance  frequency,  effective 
pavement  mass,  and  damping  constant  that  were  obtained  from  the  experi¬ 
mental  dynamic  load  deflection  curves  for  rigid  and  flexible  pavements 
axe  given  in  the  following  tabulation. 


Pavement 

Type 

Critical 

Frequency 

Hz 

Resonance 

Frequency 

Hz 

Damping 

Constant 

kip-sec/ 

in. 

Effective 
Pavement  Mass 
kip-sec^/ 
in. 

Vibrator 

Mass 

kip-sec^y 

in. 

Rigid 

concrete 

pavement 

1U.5 

10 

35 

1.0 

0.0U2 

Flexible  AC 
pavement 

15-5 

8 

8 

0.33 

0.0U2 

It.  1.1*  EFFECT  OF  STATIC  LOADS 

Some  experimental  work  has  been  done  at  WES  toward  determining 
the  dependence  of  the  dynamic  load-deflection  curves  and  the  dynamic 
stiffness  of  a  pavement -vibrator  system  on  the  static  load  applied  to 
the  pavement.  The  static  load  was  applied  to  an  AC  pavement  by  a  hy¬ 
draulic  system  that  could  apply  a  static  load  of  5-50  kips.  The  vibra¬ 
tor  was  run  at  15  Hz  for  all  of  the  tests.  Dynamic  load-deflection 
curves  were  measured  for  a  series  of  static  loads  that  was  started  at 
5  kips  and  incremented  by  adding  5-kip  lead  weights  to  the  vibrator 
system  until  a  maximum  static  load  of  50  kips  was  attained.  The  maxi¬ 
mum  value  of  the  dynamic  load  for  the  case  of  the  5-kip  static  load  was 
U  kips,  while  the  maximum  dynamic  load  for  the  case  of  the  50-kip  static 
load  was  only  10  kips.  Therefore,  for  these  experimental  tests  the  max¬ 
imum  value  of  the  dynamic  load  was  considerably  less  than  the  value  of 
the  static  load  for  which  dynamic  load-deflection  curves  were  measured. 
Figure  16  gives  the  measured  dynamic  load-deflection  curves  for  several 
values  of  the  static  load,  while  Figures  IT  and  18  show  the  dependence 
of  the  dynamic  stiffness  on  the  static  and  dynamic  loads,  respectively. 
The  experimental  data  show  that  the  dynamic  stiffness  is  greatly  depen¬ 
dent  on  the  static  and  dynamic  loads,  and  therefore  the  dynamic 
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stiffness  is  as  much  a  property  of  the  dynamic  characteristics  of  a  vi¬ 
brator  as  it  is  of  the  strength  of  the  pavement.  The  experimental  data 
in  Figure  1 6  show  that  the  dynamic  stiffness  measured  by  a  vibrator 
tends  to  increase  with  the  magnitude  of  the  applied  static  load.  The 
results  in  Figure  l6  also  show  that  the  measured  dynamic  load-deflection 
curves  tend  to  become  more  linear  with  increasing  static  load.  The 
latter  experimental  conclusion  can  be  presented  only  tentatively  because 
the  range  of  the  variation  of  the  dynamic  load  for  the  experimental 
tests  done  with  the  50-kip  _uad  was  only  0-10  kips.  The  range 

of  the  dynamic  load  should  have  been  0-40  kips  in  order  to  make  a  proper 
comparison  with  the  experimental  tests  done  with  the  5-kip  static  load 
for  which  the  dynamic  load  was  varied  from  0-4  kips. 

Figure  17  shows  that  for  the  experimental  setup  used,  the  dynamic 

stiffness  may  exhibit  a  minimum  value  for  Fs  ~  5  kips,  but  further 

tests  are  required  to  substantiate  this  conclusion.  The  results  in 

Figure  18  show  that  for  the  experimental  setup  used,  which  has  a>  <  u>c 

and  F  <  F  as  operating  conditions,  there  appear  to  be  local  minima 

s  sc 

in  the  dynamic  stiffness  as  a  function  of  F^  in  the  range  F^  ~  3-4 
kips  for  the  WES  50-kip  variable  static  load  vibrator  operating  on  AC. 

4.1.5  EFFECT  OF  BASEPLATE  RADIUS 

Experimental  work  has  been  done  at  WES  to  determine  the  dependence 
of  the  measured  values  of  the  dynamic  stiffness  on  the  size  of  the  vi¬ 
brator  baseplate.  This  experimental  work  was  done  in  part  to  confirm 
the  results  in  Section  3.7*8  regarding  the  dependence  of  the  theoreti¬ 
cally  calculated  dynamic  stiffness  on  the  baseplate  radius  of  the 
vibrator.  The  experimental  work  was  also  done  to  establish  the  exis¬ 
tence  of  the  critical  radii  (Section  3.7*7)  and  to  determine  their 
values  for  a  particular  pavement  location  so  that  the  formulas  for  cal¬ 
culating  the  layer  thicknesses  (Section  3*7*7)  and  shear  modulus  of  each 
pavement  layer  (Section  3*7*9)  could  be  verified. 

Experimental  measurements  of  the  dynamic  load-deflection  curves 
were  made  for  baseplates  having  diameters  of  3,  4,  5,  6,  7,  8,  9,  10, 

12,  15,  and  18  in.  All  of  the  baseplate  sizes  were  used  for  the  tests 
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on  rigid  pavements;  but  for  the  tests  on  AC  pavements,  the  3-,  b~,  and 
5-in.-diam  plates  could  not  be  used  because  the  pavement  failed  by 
shearing  under  the  l6-kip  static  load  of  the  vibrator.  Figure  19  gives 
the  experimental  dynamic  load-deflection  curves  for  AC  while  Figure  20 
gives  the  curves  for  the  case  of  a  rigid  pavement.  For  the  rigid  pave¬ 
ment,  the  load-deflection  curves  for  the  5-in.-diam  and  larger  baseplate 
sizes  are  very  similar,  so  that  only  the  load-deflection  curves  for  the 
5-  and  l8-in.-diam  baseplates  are  shown  in  Figure  20.  Figure  21  gives 
the  experimental  values  for  the  dynamic  stiffness  of  rigid  and  AC  pave¬ 
ments  as  a  function  of  the  baseplate  radius.  Figure  21  shows  that  an 
abrupt  change  in  the  slope  of  the  curve  giving  dynamic  stiffness  versus 
baseplate  radius  occurs  at  a  radius  of  2.5  in.  for  the  case  of  the  rigid 
pavement,  and  at  a  radius  of  6.0  in.  for  the  AC  pavement.  The  rigid 
pavement  is  known  to  have  a  thickness  of  15  in.,  while  the  AC  pavement 
is  known  to  rest  on  a  subgrade  the  top  of  which  is  at  a  depth  of  from 
28  in.  to  U2  in. 

k.1.6  COMPARISON  OF  DYNAMIC  STIFFNESS 
VALUES  MEASURED  BY  THE  DYNAFLECT 
AND  THE  WES  16-KIP  VIBRATOR 

A  comparison  of  the  dynamic  stiffness  values  determined  by  the 
WES  l6-kip  vibrator  and  by  the  Dynaflect  vibrator  on  concrete  pavements 
at  the  same  locations  is  given  in  Table  2.  Table  2  shows  that  the  stiff¬ 
ness  values  determined  by  the  Dynaflect  are  considerably  higher  than 
those  obtained  at  the  same  locations  using  the  WES  l6-kip  vibrator,  even 
though  the  radius  of  the  WES  vibrator  is  about  an  order  of  magnitude 
larger  than  the  equivalent  radius  of  the  Dynaflect  operating  on  a  con¬ 
crete  pavement  (anomalous  behavior).  The  WES  l6-kip  vibrator  has  a 
baseplate  area  of  2U5  sq  in.  The  Dynaflect  vibrator  makes  contact  with 
a  pavement  through  two  steel  wheels,  each  having  a  width  of  h  in.  and  a 
diameter  of  l6  in.  ,  and  therefore  the  contact  area  with  the  pavement 
varies  with  the  degree  of  softness  of  the  pavement,  being  considerably 
larger  for  soft  asphaltic  pavements  than  for  hard  concrete  pavements. 

The  frequency  of  operation  of  the  Dynaflect  is  8  Hz  while  that  of  the 
WES  l6-kip  vibrator  is  l6  Hz.  The  maximum  dynamic  load  of  the  Dynaflect 
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Figure  2^.  Experimental  values  of  the  dynamic  stiffness  versus  the  vibrator  baseplate  radius 
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is  0.5  kips  and  the  static  load  is  1.5  kips,  while  the  maximum  dynamic 
load  of  the  WES  l6-kip  vibrator  is  15  kips  with  a  static  load  of  16  kips. 
The  comparatively  small  size  of  the  Dynaflect  is  thought  to  be  responsi¬ 
ble  for  the  difference  between  values  of  the  dynamic  stiffness  measured 
by  the  two  vibrators  at  the  same  pavement  location.  Volume  I  of  this 
report"^  gives  a  complete  description  of  these  two  vibrators. 

4.1.7  COMPARISON  OF  DYNAMIC  STIFF¬ 
NESS  VALUES  MEASURED  BY  ROAD 
RATER  AND  WES  16-KIP  VIBRATOR 

For  the  same  locations,  dynamic  stiffness  measurements  were  made 
using  the  Road  Rater  and  the  WES  l6-kip  vibrator.  Measurements  were 
made  on  both  concrete  and  AC  pavements ,  and  the  results  are  tabulated 
in  Table  3.  As  with  the  Dynaflect,  it  was  found  that  the  Road  Rater 
generally  gives  higher  values  of  dynamic  stiffness  than  does  the  WES 
l6-kip  machine  at  the  same  locations.  This  is  true  even  though  the  WES 
l6-kip  vibrator  has  a  baseplate  radius  about  2.5  times  larger  than  the 
baseplate  radius  of  the  Road  Rater  (anomalous  behavior).  The  Road  Rater 
and  the  WES  l6-kip  vibrator  have  different  mechanical  characteristics, 
and  this  difference  is  probably  responsible  for  the  difference  in  the 
measured  values  of  the  dynamic  stiffness  at  the  same  pavement  location. 
The  Road  Rater  is  generally  operated  in  the  frequency  range  of  10-40  Hz 
and  has  a  maximum  dynamic  load  in  the  range  of  1.0  kip  and  a  static  load 
in  the  range  of  3.0  kips,  values  which  are  considerably  smaller  than  the 
corresponding  values  for  the  WES  l6-kip  vibrator.  The  Road  Rater  is 

15 

described  in  more  detail  in  Volume  I  of  this  report. 

4.1.8  COMPARISON  OF  DYNAMIC  STIFF¬ 
NESS  VALUES  MEASURED  BY  CERF 
VIBRATOR  AND  WES  16-KIP  VIBRATOR 

Table  4  gives  the  values  of  the  dynamic  stiffness  measured  at  the 
same  locations  by  the  CERF  vibrator  and  the  WES  l6-kip  vibrator.  ’ 

The  CERF  vibrator  measures  lower  values  of  dynamic  stiffness  than  does 
the  WES  vibrator.  This  behavior  of  the  dynamic  stiffness  is  in  accor¬ 
dance  with  the  larger  baseplate  contact  radius  of  the  WES  vibrator 
(normal  behavior).  The  CERF  vibrator  is  more  comparable  in  size  with 
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Table  3 

Comparison  of  Dynamic  Stiffness  Values  Measured  by  the 
Road  Rater  and  by  the  WES  l6-Kip  Vibrator 
(Flexible  Pavements) 


WES  Soil 
Test 

Stabilization 

Sections 

National  Aviation  Facilities 
Experimental  Center 

WES  l6-Kip  Vibrator  Road  Rater 

WES  l6-Kip  Vibrator 

Road  Rater 

kips/in. 

kips/in. 

kips/in. 

kips/in. 

680 

1010 

3120 

1*878 

1260 

700 

653 

2920 

1*000 

880 

653 

lll*0 

1220 

760 

1538 

21*00 

271*0 

720 

531 

2960 

5000 

2920 

381*6 

21*60 

3077 

2200 

1*31*7 

980 

960 

720 

1265 

980 

900 

800 

1123 

980 

900 

1*00 

266 

880 

862 

520 

632 

880 

833 

700 

81*0 

880 

1087 

1560 

1176 

730 

909 

I960 

2000 

900 

1087 

31*1+0 

3125 

900 

1020 

— 

— 

ON 

O 

O 

633 

— 

— 

600 

O 

GO 

LT\ 

— 

_ 

600 

667 
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Table  4 

Comparison  of  Typical  Dynamic  Stiffness  Values  as  Measured 


by  the  CERF  Vibrator  and  by  the  WES  l6-Kip  Vibrator 


Webb  Air  Force 

Base 

Kelly  Air  Force 

Base 

(Flexible  Pavement) 

(Concrete  Pavement) 

WES  l6-Kip  Vibrator 

CERF  Vibrator 

WES  l6-Kip  Vibrator  CERF  Vibrator 

kips/in. 

kips/in. 

kips/ in. 

kips/in. 

800 

680 

4400 

2725 

940 

861 

5520 

3862 

880 

783 

5440 

3026 

760 

576 

4640 

2085 

780 

696 

5840 

2627 

,  660 

587 

5440 

2875 

760 

728 

4o4o 

1957 

740 

674 

5520 

2358 

880 

779 

9200 

4895 

600 

553 

4080 

2890 

780 

746 

4920 

2174 

760 

645 

4920 

3650 

820 

657 

5600 

4000 

780 

708 

4360 

3497 

780 

698 

4560 

3129 

640 

611 

2080 

1843 

760 

636 

4800 

3109 

720 

723 

3880 

2909 

660 

664 

l64o  ■ 

1583 

960 

803 

4440 

2174 

1040 

775 

4520 

3030 

700 

469 

4o4o 

3576 

600  . 

503 

— 

— 

720 

478 

— 

— 

660 

452 
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the  WES  l6-kip  vibrator  than  is  either  the  Dynaflect  or  the  Road  Rater. 
The  CERF  vibrator  is  operated  in  the  frequency  range  of  5-60  Hz  and  has 
a  maximum  dynamic  load  of  8  kips  and  a  static  load  of  9  kips.  The  CERF 
vibrator  is  described  in  more  detail  in  Volume  I  of  this  report. ^ 

k.2  COMPARISON  OF  THEORETICAL  AND 
EXPERIMENTAL  RESULTS  FOR  THE 
DYNAMIC  STIFFNESS 

It. 2.1  DEGREE  OF  NONLINEARITY 

The  experimental  dynamic  load-deflection  curves  presented  in 
Figures  11  and  12  show  that  a  flexible  pavement  having  a  small  dynamic 
stiffness  is  more  nonlinear  than  a  rigid  pavement  with  a  large  dynamic 
stiffness.  This  experimental  fact  is  explained  in  the  selected  model  by 
the  first  nonlinear  term  in  Equation  3.22,  which  depends  on  the  dynamic 
stiffness  in  the  manner  S^  ,  so  that  the  theory  predicts  that  stiff 
pavements  can  be  expected  to  be  more  linear  than  flexible  pavements. 

It  has  been  shown  theoretically  that  for  the  case  when  w  <  uc 
and  F  <  F  ,  for  which  a  >  0  in  Equations  3.25,  3.31,  and  3.35, 

S  SC  _L 

the  dynamic  load-deflection  curves  have  a  shape  such  that  the  dynamic 
stiffness  is  a  decreasing  function  of  F^  in  the  range  of  low  values 
of  Fp  .  For  large  values  of  ,  the  dynamic  stiffness  may  continue 
to  be  a  decreasing  function  of  FD  or  may  obtain  a  minimum  value  and 
then  begin  to  increase  with  further  increase  of  the  values  of  F^  .  A 
minimum  value  for  the  dynamic  stiffness  may  or  may  not  exist  depending 
on  the  relative  values  of  the  nonlinear  parameters  b  and  e  .  Fig¬ 
ures  11  and  13  give  cases  for  which  >  0  and  for  which  the  dynamic 
stiffness  is  a  monotonic  decreasing  function  of  dynamic  load,  while 
Figures  l^t  and  18  show  cases  for  which  >  0  and  for  which  a  minimum 
occurs  in  the  value  of  the  dynamic  stiffness  at  a  given  value  of 
Fp  *  k.O  kips.  For  the  WES  l6-kip  vibrator  operated  at  a  frequency  of 
15  Hz,  it  is  found  that  the  flexible  AC  pavements  generally  have  dynamic 
load  deflection  curves  for  which  >  0  . 

It  has  also  been  shown  theoretically  that  for  the  case  when 

u  >  oj  and  F  >  F  ,  for  which  a  <  0  in  Equations  3.25,  3.31,  and 
c  s  sc  1 
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3.35.  the  dynamic  stiffness  is  an  increasing  function  of  the  dynamic 
load  in  the  region  of  low  dynamic  load.  The  hO-Hz  curve  in  Figure  lh 
may  indicate  this  behavior  with  <  0  .  Further  experimental  veri¬ 
fication  is  necessary. 

One  of  the  conclusions  of  Equation  3.31  is  that  the  dynamic  load- 
deflection  curves  measured  by  a  vibrator  will  become  essentially  linear 
in  the  region  of  small  Fp  if  the  vibrator  is  operated  at  a  character¬ 
istic  frequency  of  the  vibrator-pavement  system  which  is  called  the 
critical  frequency.  The  existence  of  the  critical  frequency  is  sug¬ 
gested  by  the  experimental  results  in  Figures  14  and  15;  however,  fur¬ 
ther  experimental  results  on  this  question  are  needed.  Figures  16  and 
17  show  conclusively  that  the  dynamic  load-deflection  curves  and  the 
dynamic  stiffness  depend  on  the  static  load  present  during  measurements. 
Although  these  curves  suggest  that  the  dynamic  stiffness  attains  a  min¬ 
imum  value  for  a  certain  value  of  the  static  load,  it  cannot  be  said 
that  this  result  has  been  definitely  verified  experimentally  and  further 
tests  are  necessary. 

h.2,2  DETERMINATION  OF  THE  PARAM¬ 
ETERS  kgg  ,  b  ,  and  e 

Only  a  limited  amount  of  numerical  analysis  has  been  done  toward 
extracting  the  parameters  kgQ  ,  b  ,  e  ,  £Q  ,  »  and  from  the 

aforementioned  experimental  data.  These  parameters  depend  on  the  radius 
of  the  vibrator  baseplate  and  on  the  elastic  parameters  v  and  G  of 
the  pavement-soil  system.  Therefore,  the  measured  values  of  these  param¬ 
eters  will  depend  on  the  type  of  vibrator  used  to  make  the  measurements 
as  well  as  on  the  elastic  parameters  of  the  pavement-soil  system. 

The  parameters  k^  ,  b  ,  and  e  can  be  obtained  by  two  methods: 

a.  If  the  dynamic  load-deflection  curves  are  sufficiently  curved 
so  as  to  allow  the  direct  application  of  a  curve-fit  program 
as  in  the  case  of  asphalt  and  weak  concrete  pavements, 
the  parameters  kQQ  ,  b  ,  and  e  can  be  obtained  diiectly 
from  the  load-deflection  curves.  A  least-square  polynomial 
fit  program  was  developed  which  fits  the  dynamic  load- 
deflection  curves  with  an  odd  polynomial  containing  fifth-, 
third-,  and  first-order  terms.  This  fitting  procedure  immedi¬ 
ately  determines  ,  -sa.^  ,  and  that  appear  in 
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Equation  3.22.  The  value  of  kg  is  then  obtained  from  Sq 
by  using  Equation  3.21  and  appropriate  values  for  the  pave¬ 
ment  mass  and  damping  constant.  Placing  the  values  of  , 
<*2  »  Sq  ,  and  kg  into  Equations  3.25  and  3.26  immediately 
gives  the  values  of  b  and  e  .  The  values  of  kgg  and  x 
are  then  obtained  from  a  simultaneous  solution  of  Equations 
3.1*  and  3.15. 

b.  The  dynamic  load-deflection  curves  for  very  rigid  pavements 
are  essentially  straight  lines.  The  dynamic  load-deflection 
curves  are  measured  for  a  series  of  static  loads  by  using  a 
vibrator  with  a  variable  static  force.  The  dynamic  stiffness 
is  determined  from  the  slope  of  the  dynamic  load-deflection 
curves.  Because  the  slope  is  essentially  independent  of  the 
dynamic  load,  it  follows  from  Equation  3.23  that 
S  ^  Sg(uj,F  )  ,  so  that  the  measured  dynamic  stiffness  value 
immediately  determines  Sg(w,Fs)  .  The  function  kg(Fs)  is 
then  determined  from  Sg(u>,Fs)  by  using  Equation  3.21.  The 
function  kg(Fg)  is  then  fitted  to  the  expression  in  Equa¬ 
tion  3.15  after  x^  has  been  expressed  in  terms  of  Fs  by 
using  Equation  3.2.  This  fitting  procedure  determines  k  , 
b  ,  and  e  .  00 

Method  a  can  be  used  for  those  pavements  which  have  dynamic  load- 
deflection  curves  that  are  sufficiently  curved  to  be  able  to  obtain 
al  ’  “p  »  and  directly  from  a  least-squares  program.  This  is  not 

generally  the  case  for  very  strong  PCC  pavements.  For  very  rigid  pave¬ 
ments,  method  b  should  be  used  because  the  dynamic  load-deflection 
curves  are  essentially  straight  lines.  Computer  programs  have  been 
developed  which  calculate  the  parameters  kQQ  ,  b  ,  and  e  directly 
from  the  measured  dynamic  load  deflection  curves  using  method  a.  No  nu¬ 
merical  work  has  been  done  applying  method  b  to  the  rigid  pavement. 
Figures  22,  23,  and  2h  give  the  values  of  kQQ  ,  b  ,  and  e  ,  respec¬ 
tively,  that  were  determined  by  method  a  from  the  measured  dynamic  load- 
deflection  curves  for  a  series  of  baseplate  radii.  The  existence  of  the 
first  critical  radius  can  be  established  from  the  abrupt  change  in  slope 
of  the  function  k^Ca)  for  a  =  .  Method  a  was  applied  to  the  dy¬ 

namic  load-deflection  curves  of  a  rigid  pavement  and  a  flexible  pavement. 
The  considerable  scatter  in  the  results  for  b  and  e  for  the  rigid 
and  flexible  pavements  probably  reflects  the  difficulty  of  determining 
the  curvature  of  the  dynamic  load-deflection  curves  which  may  contain 
some  experimental  error.  The  large  scatter  in  the  experimentally 
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Figure  22.  Experimentally  derived  values  of  the  parameter 
kQQ  versus  baseplate  radius 


lAMPlATC  RADIUS,  •  A 

Figure  23.  Experimentally  derived  values  of  the  parameter 
b  versus  baseplate  radius 
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BASEPLATE  RADIUS,  IN. 

b>;  FLEXIBLE  PAVEMENT 

Figure  2b  Experimentally  derived  values  of  the  parameter 
e  versus  baseplate  radius 
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derived  values  of  the  parameters  b  and  e  probably  means  that  only 
the  experimentally  derived  values  of  kgg  for  a  series  of  baseplate 
sizes  will  be  useful  for  calculating  pavement  structure. 

The  values  of  cm  be  calculated  from  Equations  3.103-3  107 
provided  the  Poisson's  ratio  of  each  layer  of  the  pavement  system  is 
known  and  provided  the  values  of  the  critical  radii  are  known.  The 
values  of  and  can  be  obtained  in  terms  of  the  values  of  b 

and  e  given  in  Sections  3. 7. 2-3. 7*5*  "  \dded  the  shear  modulus  and 
layer  thickness  are  known  for  each  laye*.  contained  within  the  static 
finite  depth  of  influence.  The  experimental  values  of  ,  and 

appear  in  Figures  25,  26,  and  27,  respectively.  The  results  in 
Figure  25  show  that  is  larger  for  rigid  than  for  flexible  pave¬ 
ments.  The  large  scatter  in  the  values  of  ln  and  are  probably 

due  to  experimental  variations  in  the  curvature  of  the  dynamic  load- 
deflection  curves.  Typical  values  of  kQQ  ,  b  ,  e  ,  ,  SL^  »  and 

for  the  WES  l6-kip  vibrator  with  an  l8-in.-diam  baseplate  are  given 
in  Table  5. 

4.2.3  DETERMINATION  OF  THE  SHEAR 
MODULI  AND  THICKNESSES  OF 
PAVEMENT  LAYERS 

The  thickness  of  each  pavement  layer  can  be  determined  using 
Equations  3.108-3.110,  and  the  shear  modulus  of  each  layer  can  be  de¬ 
termined  by  using  Equations  3.117-3.123.  These  equations  require  Pois¬ 
son's  ratio  of  each  pavement  layer,  the  critical  radii,  and  the  function 
kgg(a)  as  input  variables.  The  critical  radii  are  determined  by  measur¬ 
ing  the  dynamic  load-deflection  curves  for  a  series  of  baseplate  sizes 
and  noting  the  particular  values  of  the  baseplate  radius  at  which  abrupt 
changes  occur  in  the  slope  of  the  curve  showing  dynamic  stiffness  versus 
baseplate  radius.  The  critical  radii  can  also  be  obtained  by  noting  the 
values  of  the  baseplate  radius  at  which  abrupt  changes  in  the  slope  of 
kgg(a)  occur.  Table  6  gives  the  values  of  critical  radii,  layer  thick¬ 
nesses,  and  shear  moduli  corresponding  to  the  experimental  data  shown  in 
Figure  21.  Only  the  first  critical  radius  was  measured  experimentally 
for  the  rigid  pavement,  and  only  the  second  critical  radius  was  observed 
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Figure  25.  Experimentally  derived  values  of  the  parameter 
versus  vibrator  baseplate  radius 


Figure  26.  Experimentally  derived  values  of  the  parameter 
versus  baseplate  radius 


WES  TEST  SECTIONS 


t  »K^D  PAVEMENT 
A  flllllLl  PAVEMENT 


BASEPLATE  RADIUS.  IN. 


Figure  27.  Experimentally  derived  values  of  the  parameters 
versus  baseplate  radius 


Table  5 

Approximate  Numerical  Values  of  the  Parameters 


Appearing  in  the  Nonlinear  Pavement  Model 


WES  l6-Kip  Vibrator’ 


Parameter 


Concrete  Pavement  Asphaltic  Concrete  Pavement 


8.0  x 

10°  lb/in. 

2.0  x 

10  lb/ in. 

XT' 

O 

X 

109  lb/in.3 

-1.3 

x  109  lb/in. 

X 

o 

aj 

10ll+  lb/in.5 

1*. 0  x 

1012  lb/in. 

200  in. 

72  in. 

1.0  x 

105  in."1 

2.0  x 

1014  in.-1 

3.0  x 

1010  in.'3 

-3.5  x 

10T  in."3 

Table  6 

Experiments]  Values  of  the  Parameters  Describing 
the  Subsurface  Structure  of  Pavements 
at  the  WES  Test  Area 


Parameter 

Ritfid  Pavement 

Flexible  Pavement 

First  critical  radius,  a  , 

cl 

2.5  in. 

Not  measured 

Second  critical  radius,  a  _ 

c2 

Not  measured 

6.0  in. 

Thickness  of  pavement,  h^ 

15.0  in. 

33.0  in. 

Poisson's  ratio, 

0.25 

0.35 

Shear  modulus , 

2.U  x  10^  psi 

3.7  x  10U  psi 

Compression  modulus,  E^ 

6. 0  x  10^  psi 

1.0  x  10^  psi 

Volume  factor  of  the  frustum 
of  cone  of  stress,  ¥ 

3.8 

3.8 

Ratio  of  radii  of  frustum 
of  cone  of  stress,  k 

2.8 

2.8 

for  the  flexible  pavement.  The  determination  of  the  thickness  of  a  thin 
upper  layer  of  a  rigid  or  flexible  pavement  requires  the  use  of  very 
small  baseplate  sizes.  The  shear  failure  of  the  flexible  pavement 
prevented  any  tests  using  the  WES-16  kip  vibrator  for  a  baseplate  whose 
diameter  was  smaller  than  3  in. ,  and  this  prevented  the  measurement  of 
the  first  critical  radius  for  the  flexible  pavement  at  the  WES  test 
section. 

The  experimental  data  presented  in  Section  U.1.5  generally  tend 
to  confirm  the  existence  of  the  critical  radii  that  were  introduced  in 
Section  3. 7*7.  For  the  case  of  the  rigid  pavement  shown  in  Figure  21, 
the  abrupt  change  in  the  slope  of  the  dynamic  stiffness  function  for  a 
baseplate  radius  of  2.5  in.  is  interpreted  to  mean  that  the  value  of  the 
critical,  radius  is  2.5  in.  For  this  value  of  the  baseplate  radius,  the 
zero1'  1  term  of  the  finite  depth  of  influence,  ,  as  defined  by  Equa¬ 
tions  3.103  and  3.108,  is  equal  to  the  depth  of  the  first  interface 
which  in  this  case  is  15  in.  The  condition  =  15  in.  Combined 
with  Equations  3*52,  3.103,  and  3.108  give  the  values  of  V  and  k  that 
appear  in  Table  6.  The  shear  modulus  of  the  upper  pavement  layer  is 
calculated  from  the  slope  of  the  experimentally  derived  function  kQ0(a) 
that  appears  in  Figure  22  by  using  Equation  3.117. 

h.2.h  ANOMALOUS  VALUES  OF  DYNAMIC 
STIFFNESS 

The  difference  in  the  values  of  the  dynamic  stiffness  measured  by 
the  Dynaflect,  Road  Rater,  CERF  vibrator,  and  the  WES  l6-kip  vibrator  at 
the  same  pavement  location  is  thought  to  be  due  to  the  aifferent  base¬ 
plate  size,  static  load,  and  dynamic  load  used  by  these  vibrators.  The 
difference  in  the  values  of  the  measured  dynamic  stiffness  can  be  con¬ 
siderable  as  has  been  described  in  Sections  1*.  1.6-4. 1.8  and  in  Tables  2- 
U.  The  nonlinear  mechanical  model  of  the  dynamic  response  of  pavements 
to  surface  loads  that  was  developed  in  Section  3  implies  that  the  mea¬ 
sured  values  of  dynamic  stiffness  should  depend  on  the  baseplate  radius, 
static  load,  and  dynamic  load  of  the  vibrator  that  is  used  to  make  these 
measurements.  This  nonlinear  mechanicai.  model  should  be  sufficient  to 
correlate  the  different  values  of  dynamic  stiffness  measured  by 
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different  sized  vibrators  at  the  same  pavement  location.  However,  no 
numerical  studies  of  this  correlation  have  been  performed,  and  only 
qualitative  results  will  be  described. 

Dynamic  stiffness  measurements  were  made  at  the  same  location 
using  different  vibrators  in  order  to  determine  further  the  dependence 
of  the  measured  values  of  the  dynamic  stiffness  on  the  type  of  vibrator 
used  to  make  the  measurements.  These  measurements  produced  the  unusual 
result  that  under  some  circumstances  it  is  possible  at  the  same  pavement 
site  to  measure  a  larger  value  of  dynamic  stiffness  with  a  small  base¬ 
plate  than  with  a  large  baseplate  (Sections  4.1.8-4.1.10  and  Tables  4-6). 
The  nonlinear  pavement  model  that  was  developed  in  Sections  3. 5-3. 8  at¬ 
tempts  to  explain  these  anomalous  experimental  results  for  the  dynamic 
stiffness  measurements. 

For  flexible  pavements,  a  small  part  of  the  discrepancy  between 
the  values  of  dynamic  stiffness  measured  by  the  WES  l6-kip  vibrator  and 
those  measured  by  the  Dynaflect,  Road  Rater,  and  CERF  vibrators  at  the 
same  location  is  due  to  the  fact  that  the  values  of  dynamic  stiffness 
presented  for  the  WES  l6-kip  vibrator  in  Volume  I  of  this  report1'*  were 
determined  from  the  slope  of  the  dynamic  load-deflection  curves  in  the 
region  of  large  dynamic  load  (*“12  kips).  These  values  of  dynamic  stiff¬ 
ness  are  considerably  lower  than  the  values  that  would  be  obtained  from 
the  region  of  lower  dynamic  load  (Figure  11).  Better  agreement  with  the 
Dynaflect,  Road  Rater,  and  CERF  vibrators  would  be  obtained  for  flexible 
pavements  if  the  values  of  dynamic  stiffness  measured  by  the  WES  l6-kip 
vibrator  were  obtained  from  that  portion  of  the  dynamic  load-deflection 
curves  where  the  dynamic  load  is  comparable  to  the  lower  value  produced 
by  the  other  vibrators  in  question.  Another  contribution  to  the  dis¬ 
crepancy  between  the  values  of  dynamic  stiffness  measured  by  different 
vibrators  at  the  same  locations  is  due  to  the  difference  in  operating 
frequency  as  indicated  in  Table  1. 

Actually,  the  dependence  of  the  measured  values  of  the  dynamic 
stiffness  on  the  magnitude  of  the  dynamic  load  of  the  vibrator  plays 
only  a  part  in  the  explanation  of  the  anomalous  dynamic  stiffness  mea¬ 
surements.  This  explanation  of  the  anomalous  values  of  the  dynamic 
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stiffness  measurements  does  not  apply  to  strong  rigid  pavements  whose 
dynamic  load-deflection  curves  are  essentially  linear  (Figure  12).  In 
this  case  the  dynamic  stiffness  is  essentially  independent  of  the  mag¬ 
nitude  of  the  dynamic  load,  and  the  explanation  of  the  higher  values  of 
dynamic  stiffness  measured  by  a  vibrator  with  a  small  baseplate  contact 
radius  (about  1  in.  in  the  case  of  Dynaflect)  compared  to  those  measured 
using  a  vibrator  with  a  large  baseplate  contact  radius  (about  9  in.  in 
the  ease  of  the  WES  l6-kip  vibrator)  depends  on  the  two  nonlinear  ef¬ 
fects  shown  in  Figure  4a  and  described  in  Section  3.8.2: 

a.  The  spring  constant  k_  is  a  function  of  F^  which  has  a 
local  minimum  value  for  a  specific  value  of  sFs  . 

b.  The  parameter  k^  is  not  a  linear  function  of  baseplate 
contact  radius  for  the  case  of  a  layered  pavement  structure 
which  has  a  finite  depth  of  influence,  but  is  a  very  slowly 
increasing  function  of  baseplate  radius  when  the  baseplate 
radius  is  larger  than  the  first  critical  radius. 

4.2.5  CALCULATION  OF  AVERAGED 
PARAMETERS  FOR  A  NON- 
HOMOGENEOUS  PAVEMENT 

The  static  and  dynamic  stress  and  strain  distribution  in  the 
pavement  beneath  the  vibrator  baseplate  is  three-dimendional  in  nature. 

The  average  values  that  are  calculated  in  Equation  3.66  are  based  on  an 
assumed  one-dimensional  distribution  of  strain  in  the  pavement.  Some 
experimental  evidence  on  the  relatively  large  surface  deflection  of  pave¬ 
ments  associated  with  an  applied  static  load  suggests  that  the  average 
values  that  appears  in  Equation  3.66  should  more  properly  be  volume  aver¬ 
ages.  To  account  for  the  volume  averages,  it  is  a  simple  matter  to  in¬ 
sert  volume  ratio  factors  for  each  layer  thickness  that  appears  in  Equa¬ 
tion  3.66.  Each  pavement  layer  is  associated  with  a  strain  distribution 
which  can  be  represented  by  a  frustum  of  a  cone.  The  ratios  of  the 
volume  of  the  frustum  of  the  cone  in  each  pavement  layer  to  the  volume 
of  the  complete  frustum  of  the  cone  of  strain  in  the  pavement  and  sub¬ 
grade  are  the  volume  ratio  factors  that  should  replace  the  one-dimensional 
average  that  appears  in  Equations  3.66-3.99*  The  volume  average  tends  to 
increase  the  effect  of  the  subgrade  as  compared  to  that  of  the  upper 
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5.  SUMMARY,  CONCLUSIONS,  AND  RECOMMENDATIONS 


5.1  SUMMARY 

A  nonlinear  mechanical  model  describing  the  dynamic  properties  of 
a  pavement -vibrator  system  has  been  developed  which  describes  the  non¬ 
linear  dynamic  response  of  a  pavement  to  a  sinusoidal  loading  applied  to 
the  pavement  surface.  Theoretical  expressions  are  developed  for  the  dy¬ 
namic  stiffness  of  a  pavement  measured  by  a  mechanical  vibrator  which 
are  expressed  in  terms  of  the  static  load,  dynamic  load,  and  vibrator 
baseplate  size  and  in  terms  of  the  linear  and  nonlinear  pavement  param¬ 
eters.  The  nonlinear  mechanical  model  gives  an  analytical  correlation 
among  the  values  of  the  dynamic  stiffness  measured  by  different  vibra¬ 
tors  at  the  same  pavement  location.  The  nonlinear  parameters  that 
characterize  the  dynamic  pavement-restoring  force  are  calculated  in 
terms  of  the  layered  structure  of  the  pavement-soil  system  and  in  terms 
of  the  finite  depth  of  influence  of  the  static  load  of  the  vibrator. 

The  basic  nonlinearity  of  the  pavement  system  is  used  to  determine  the 
shear  modulus  and  thickness  of  each  pavement  layer  from  the  measured 
values  of  the  dynamic  stiffness  for  a  series  of  baseplate  sizes.  Ex¬ 
perimental  tests  were  done  to  determine  the  validity  of  the  theoretical 
pavement  response  model.  The  experimental  and  theoretical  results  are 
in  good  agreement. 

5.2  CONCLUSIONS 

The  nonlinear  mechanical  model  developed  in  this  report  gives  the 
following  conclusions: 

a.  Third-  and  fifth-order  nonlinear  terms  in  the  displacement 

are  required  to  describe  the  dynamic  load-deflection  response 
of  actual  pavement  systems  (Section  3.2). 

b_.  The  theoretical  nonlinear  oscillator  model  of  pavement  re¬ 
sponse  to  a  dynamic  loading  shows  that  stiff  pavements  have 
a  more  linear  dynamic  load-deflection  curve  than  flexible 
pavements  (Section  3.5.3). 

c_.  At  specific  critical  frequencies  the  dynamic  load-deflection 
curves  become  essentially  linear  at  low  values  of  the  dynamic 
force  (Section  3.5«**). 


d.  An  explanation  of  the  nonlinear  response  of  pavement  and  soil 
systems  requires  a  departure  from  linear  elasticity  and  a 
finite  depth  of  influence  for  a  static  load  applied  to  the 
surface  of  a  pavement  (Section  3.6). 

e.  The  linear  and  nonlinear  parameters  of  the  nonlinear  pavement 
restoring  force,  k  ,  b  ,  and  e  ,  can  be  described  in 
terms  of  the  elastic  constants  and  thicknesses  of  the  pave¬ 
ment  layers.  In  general,  b  and  e  are  discontinuous  func¬ 
tions  of  the  vibrator  baseplate  radius  and  have  the  proper¬ 
ties  b  <  0  and  e  >  0  ,  while  k  is  a  continuous  and 
increasing  function  of  the  baseplate  radius  (Section  3.7). 

f.  The  anomalous  values  of  the  dynamic  stiffness  wherein  a  vi¬ 
brator  with  a  small  baseplate  and  a  small  static  load  can 
measure  a  larger  dynamic  stiffness  value  them  a  vibrator  with 
a  large  baseplate  and  a  large  static  load  is  explained  in 
terms  of  the  penetration  of  the  finite  depth  of  influence 
through  the  successive  layers  of  the  pavement  system  and  in 
terms  of  the  dependence  of  the  dynamic  stiffness  on  the 
static  load  (Section  3.8.2). 

£.  The  measured  value  of  the  dynamic  stiffness  is  not  entirely 
a  property  of  the  pavement  because  it  also  depends  on  the 
size  of  the  dynamic  load,  static  load,  and  baseplate  radius 
of  the  vibrator  that  is  used  to  make  the  measurement. 
(Sections  3-^  and  3.5). 

h.  Critical  values  of  the  baseplate  radius  of  the  vibrator  are 
expected  to  exist  for  which  the  dominant  term  of  the  static 
finite  depth  of  influence  passes  through  the  various  inter¬ 
faces  of  the  pavement  layers.  Abrupt  changes  in  the  slope 
of  the  dynamic  stiffness  versus  baseplate  radius  "function  are 
expected  to  occur  at  the  critical  radii  (Sections  3.7.7  and 
3.7.8). 

i_.  The  variation  of  the  dynamic  stiffness  with  the  vibrator 
baseplate  radius  is  less  rapid  than  is  predicted  by  the 
linear  elastic  theory.  For  rigid  pavements  the  dynamic  stiff 
ness  is  roughly  independent  of  baseplate  radius  when  the  base 
plate  radius  is  larger  than  the  first  critical  radius  (Sec¬ 
tions  3-7.8  and  3.8.2). 

J_.  A  determination  of  the  elastic  constants  of  each  layer  in  a 
pavement  system  may  be  possible  from  dynamic  stiffness  mea¬ 
surements,  but  requires  that  the  vibrator  have  the  capability 
of  varying  the  radius  of  its  baseplate  (Sections  3.7.7,  3.7.8 
and  3.7.9). 

It.  The  thickness  of  each  pavement  Txver  can  be  calculated  in 
terms  of  the  critical  radii  (Section  3.7.7). 

A  comparison  of  the  theoretical  and  experimental  results  for  the 
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dynamic  load-deflection  curves  and  the  dynamic  stiffness  yields  the 
following  conclusions: 

fi.  Some  of  the  experimental  dynamic  load-deflection  curves  of 
pavements  are  nonlinear.  This  is  especially  true  of  more 
flexible  AC  pavements.  The  stiffer  PCC  pavements  exhibit 
dynamic  load-deflection  curves  that  are  only  slightly  non¬ 
linear  (Section  U.l). 

b.  Generally  speaking,  stiff  pavements  have  dynamic  load  deflec¬ 
tion  curves  which  are  more  linear  than  those  corresponding  to 
more  flexible  pavements  (Section  U.1.2  and  Figure  11). 

c_.  There  is  some  evidence  that  the  critical  frequency  exists  and 
that  it  has  a  value  of  about  15  Hz  for  AC  and  PCC  pavements 
(Sections  3.5.1*,  3-5.5,  and  U.1.3  and  Figures  11,  lU,  and  15) 

d.  The  value  of  the  dynamic  stiffness  of  a  given  pavement  de¬ 
pends  on  the  magnitude  of  the  dynamic  force  static  load,  and 
baseplate  size  of  the  vibrator  (Sections  3.1+,  3.5,  3.8,  U.1.2 
l+.l.U,  1.5,  Figures  16-19,  and  Tables  2-k). 

£.  The  measurement  of  dynamic  load-deflection  curves  for  a  se¬ 
ries  of  baseplate  sizes  allows  the  critical  radii  to  be  de¬ 
termined.  The  shear  modulus  and  thickness  of  each  pavement 
layer  can  be  determine,’  in  terms  of  the  critical  radii,  the 
measured  values  of  dynamic  stiffness  in  the  intervals  between 
the  critical  radii,  and  the  assumed  value  of  Poisson's  ratio 
for  each  pavement  layer  (Sections  3.7-7,  3.7.9,  and  U.1.7). 

5.3  RECOMMENDATIONS 

The  positive  conclusions  of  this  report  regarding  the  usefulness 
of  a  mechanical  vibrator  for  determining  the  strength  of  pavements  and 
their  subsurface  structure  suggest  that  a  more  detailed  study  of  the 
measured  dynamic  stiffness  function  be  made. 

5.3.1  DETERMINATION  OF  SUBSURFACE  STRUCTURE 

In  this  report  only  one  case  of  a  rigid  pavement  and  one  case  of 
a  flexible  pavement  were  used  for  the  determination  of  the  shear  modulus 
and  thickness  of  each  pavement  layer  directly  from  the  values  of  the  dy¬ 
namic  load-def lection  curves  that  were  measured  at  the  pavement  surface 
using  a  mechanical  vibrator  with  a  variable  baseplate  size.  It  is  sug¬ 
gested  that  many  different  test  sites  of  known  subsurface  structure  be 
tested  with  the  purpose  of  validating  the  method  of  subsurface  structure 


106 


determination  which  employs  a  variable  baseplate  size  that  is  presented 
in  Sections  3*7*7  and  3*7*9* 


i 


i 


; 

{ 


* 


5*3*2  CORRELATION  OF  DYNAMIC  STIFF¬ 
NESS  VALUES  MEASURED  BY 
DIFFERENT  VIBRATORS 

This  report  has  determined  the  theoretical  dependence  of  the  dy¬ 
namic  stiffness  values  on  the  size  of  the  dynamic  load,  static  load,  and 
baseplate  radius  that  are  used  during  the  measurement  of  the  dynamic 
stiffness  of  a  pavement.  This  theoretical  dependence  enables  the  values 
of  the  dynamic  stiffness  measured  by  different  vibrators  to  be  compared. 
A  more  detailed  numerical  study  should  be  performed  to  validate  the  ca¬ 
pability  of  analytically  correlating  the  values  of  the  dynamic  stiffness 
measured  by  different  vibrators  at  the  same  pavement  location. 

5*3*3  SELECTION  OF  VIBRATOR 

The  mechanical  vibrators  that  are  used  to  evaluate  pavements  non- 
destructively  should  be  viewed  as  scientific  instruments  that  are  used 
to  determine  the  nature  and  structure  of  pavement  systems.  Each  pave¬ 
ment  and  its  corresppnding  operational  traffic  will  require  its  own 
specific  vibrator  characteristics  if  the  nondestructive  testing  is  to 
be  accomplished  optimally.  From  the  conclusions  of  this  report  it  fol¬ 
lows  that  optimal  operating  conditions  can  be  obtained  only  if  the  vi¬ 
brator  has  the  capability  of: 

a.  Applying  a  series  of  static  loads. 

b.  Generating  a  series  of  dynamic  loads  over  a  range  of  constant 
frequencies. 

c_.  Operating  in  a  frequency  range  which  includes  the  critical 
frequency. 

4.  Applying  a  series  of  baseplate  sizes  to  the  pavement. 

It  is  recommended  that  the  construction  of  a  new  vibrator  should  include 
the  capability  of  easy  changes  of  baseplate  size  and  magnitude  of  the 
static  and  dynamic  loads  applied  to  the  pavement  surface. 
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